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DISCLAIMER 



Although each program has been tested by its contributor, no 
warranty, express or implied, is made by the contributor or 
1620 USERS Group, as to the accuracy and functioning of the 
program and related program material, nor shall the fact of 
distribution constitute any such warranty, and no responsibility 
is assumed by the contributor or 1620 USERS Group, in con- 
nection therewith. 
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Title : 

"Solution of Homogeneous and Non-Homogeneous Simultaneous 
Linear Equations" 

Subject Classification : 

5.0 

Author; Organization : 

Heberto Pachon - Staff Mathematician 
Automatic Electric Laboratories, Inc. 

Direct Inquiries to : 

Heberto Pachon, Automatic Electric Laboratories, Inc., 400 
North Wolf Road, Northlake, Illinois; Telephone: 345-7111 ext. 587. 

Purpose/Description : 

The program solves homogeneous and non-homogeneous systems of 
linear equations; at the same time, it provides a measure of the 
accuracy of the computed solution. 

If the system has a unique solution, this solution is typed. 
If the system has an infinite number of solutions, a general solu- 
tion can be obtained. 

If the system is inconsistent, it has no solution and this is 
the answer typed. 

Mathematical Method : 

The system' s augmented matrix formed by the system 1 s coefficient 
matrix and the constant column matrix, is reduced by means of elemen- 
tary linear operations to an equivalent matrix, in which the system's 
coefficient matrix appears as a diagonal matrix with only 1 and in 
the main diagonal. If only l's appear in the main diagonal, the sys- 
tem has a unique solution; if l's and 0's appear in the main diagonal, 
the system may have an infinite number of solutions or no solution at 
all, depending on the values of the equivalent constant column matrix. 



Restrictions, Ran^e: 

The program is restricted to "Non-ill-conditioned" systems, 
and to a maximum of 17 equations with 17 unknowns. 

Storage Requirements: 

19938 positions 
Equipment Specifications: 

Memory 20K X 40K 60K K Automatic Divide: Yes X No 

Indirect Addressing: Yes X N o Other Special Features 

Required : Card system 

Additional Remarks: 

(a) Program is written in FORTRAN with FORMAT. 

(b) The program has run successfully about 20 times; these 
include solution of a maximum system (17 equations with 17 unknowns). 

(c) The program was used successfully in conjunction with the 
calculation of the gain, input impedance, and output impedance of a 
transistorized amplifier designed in our Laboratories. 

(d) If source deck is compiled to obtain Program deck, the 
FORTRAN with FORMAT system used must contain the ABSOLUTE VALUE 
SUBROUTINE. 

(e) If Program deck is used directly, it already contains the 
ABSOLUTE VALUE SUBROUTINE. 
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Deck Labeling 



Two decks are being submitted. 



Deck #1 includes a Hash total card, three comment cards, and all the 
FORTRAN statements comprising the program. When the Hash total card 
is removed from this deck, the deck left is known as the SOURCE deck. 
This deck can be compiled by the FORTRAN with FORMAT processor in 
order to obtain what from now on will be called the PROGRAM deck. 



Deck #2 includes a Hash total card, and 582 cards containing the 
machine language instructions comprising the program. 

When the Hash total card is removed from this deck, the deck left 
is known as the PROGRAM deck. 
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Program Description 

The aim of this program is to provide an approximate solution to a 
homogeneous , or non -homogeneous , system of linear equations. The 
program is restricted to a maximum system of 17 equations with 17 
unknowns, and to "non-ill-conditioned" 1 systems. The question of 
extending the program's range to systems which contain a larger 
number of equations, or to "ill-conditioned" systems, is out of 
direct consideration. This restriction is imposed here because the 
program is set up for single precision arithmetic in a computer which 
carries only eight significant digits on any arithmetic operation, 
and has an available 20K memory. With a larger core storage, the 
program could be set up for multiple precision arithmetic, and this 
is at the present time the most common way to extend the accuracy 
and the size of the system that can be handled. 

The program provides core storage for a rectangular matrix of 17 x 18 
maximum size, but it will adjust itself for systems whose augmented 
matrix is smaller than this size. 

For any N x N system (N s 17) , the program can obtain solutions for 
as many as 18-N constant vectors which carry the name "VECTOR N+l, 
VECTOR N+2, , VECTOR 18." 

When a solution has been determined, the sentence SOLUTION FOR VECTOR 
L is typed; under this heading the appropriate solution is typed. 

If the equations are independent of each other, there is a unique 
solution for the system, and the solution values are typed one per 
line. These values correspond to the unknowns X x , X a , ,X n . 

If the equations are not all independent, and together with the 
constant vector they form a consistent system, there is an infinite 
number of solutions for the system, and a message is typed to indicat 
so. This message also provides two choices; to obtain a general solu 
tion, or to continue with the next problem. 

If the equations are not all independent, and together with the con- 
stant vector they form an inconsistent system, there is no solution 
for the system, and a message is typed to indicate inconsistency. 



1 See "Determination of the Zero Value" section, for explanation 
of meaning. 
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The program is based on the theory described under Method of Cal- 
culation, and consists of reducing the system's augmented matrix 
to the row equivalent matrix described in Theorem 1. Analysis of 
the elements forming this equivalent matrix determines whether the 
system is consistent or not. 

Two problems arise at this time. First, when checking for incon- 
sistency or linear dependence, elements of the equivalent matrix in 
Theorem 1 must be compared against a computed zero value. Actual 
computations showed that a computed zero very seldom reaches the 
value zero; this is due to the round-off and truncation errors in- 
troduced by arithmetic computations. Second, when a solution has 
been reached, how accurately do the calculated unknowns satisfy 
the system under solution. 

In order to solve these problems, a tolerable zero value i s introduced 
in the program, and a rule to determine this value is explained under 
Determination of the Zero Value. This tolerable zero value worked out 
for all examples analyzed by the program, and it seems to be a good 
approximation; however, the author does not imply that this particular 
approximation is absolutely correct. 

This zero value is also a measure of the accuracy of the solution 
obtained. In other words, once the unknowns are typed out, an 
evaluation of the left side of an equation in the system will show 
that it differs from its equality constant by this zero value, at 
most. 

If the calculated unknowns had not satisfy the equations within the 
accuracy specified by the zero value, they would not have been typed. 
Instead, a message suggesting to increase the value of zero would 
have been typed . 

This accuracy self -adjustment provided for in the program and the 
program's ability to determine whether the system has a unique solu- 
tion, an infinite number of solutions, or no solutions at all, are 
features not provided for by any of the existing programs that the 
author has knowledge of. 

Another extra feature is the program's ability to generate the general 
solution to a system containing an infinite number of solutions. Such 
is the case in a homogeneous system of two equations with three un- 
knowns, for example. 
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Once the object deck has been obtained and loaded, any number of 
systems can be solved one after the other without reloading the 
Program Deck. 

Additional information about the characteristics of the program is 
contained in the sections to follow. 



Method of Computation 

The computations carried out by this program are based on the theorems 
appearing below. These theorems are fully explained and proved by 
S. Perlis 1 , or F. B. Hildebrand a , and any one interested in a deeper 
mathematical background of the material presented here should refer to 
their books. The theorems are listed here for a better understanding 
of the program. 

Theorem 1 - Every n x m rectangular matrix C can be reduced to a 
matrix D with the following properties: 

(1) For some r s 0, all elements d^j (j_ = r+l] x+2,™. . . . , n) 
are equal to zero. ' ' ' 

(2) The first non-zero element appearing in the K th row 
(1 s K £ r) is equal to 1; or d K j = 1 for some j - L 
( 1 £ L £ m) , and d k j - (j < L, L ^ 1) . The column 
L in which this d k j occurs being numbered CL. 

(3) All elements on column CL above are equal to zero, except 
the unity element in row K. 

(4) There exists exactly r columns of the type described in 3, 
or 1 £ L £ r, and Cj < C 3 <.....< C r . 

D is said to be "row equivalent to C." An example of a matrix with 
properties 1-4 above is : 



1 Sam Perlis, Theory of Matrices , (Reading, Massachusetts: Addison- 
Wesley Publishing Company, Inc., 1958), pp. 36-48. 

• F. B. Hildebrand, Methods of Applied Mathematics , (Englewood Cliffs, 
N.J.: Prentice-Hall, Inc., 1961), pp. 18-23. 
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D - 



'1 5 5 

1 2 4 

1 3 







2 10 10 

3 15 1 2 19 

4 20 3 6 1 35 

5 25 4 8 3 50 

6 30 5 10 4 62 



for which r - 3, 1 £ L £ 3, and C x , C a , and C 3 correspond to columns 
number one, three, and five, respectively. 

The steps followed in the reduction of C to D are called "elementary 
operations." 

Theorem 2 - If A and are the matrices resulting from deleting the 
extreme righthand column of the matrices C and D of Theorem 1, respec- 
tively, then Aj is row equivalent to A and A x has the properties 1-4. 

By definition, the "row rank" of a matrix C is the maximum number of 
linearly independent rows in C. With this definition, the proof of 
the following theorem is obvious . 

Theorem 3 - The row rank of a matrix C is the number "r" of non-zero 
rows appearing in its equivalent matrix D of Theorem 1. 

Now, if we consider a system of n linear equations in unknowns 
Xj , X, , , Xjj of the form 

a 11 X x + a 13 Xjj + + a in X n = k x 

a ai Xj + a aa X 3 + + a an X n = 

a n i X t a n8 X a + + a nn Xn " k n 

and let 

A - (aij) (n x n) matrix 

X = col (X 1} Xg, •, X n ) (n x 1) matrix 

K - col (k x , k 8 , , k n ) (n x 1) matrix 
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the system becomes 



AX = K 



which is the common matrix notation of a system of linear equations, 
and in which A is called the "coefficient matrix" of the system. 

Going one step further, we lee 

C = (A, K) (n x (n + 1)) matrix 

be the "augmented matrix" of the system. Then, Theorem 1 indicates 
that C can be reduced, to an equivalent matrix 

D = (A t , 1^ ) 

Using Theorem 2, we conclude that A *» A x , which indicates that AX = K 
is equivalent to A t X = 1^ . 

Let us assume now that once D is found, by Theorem 3 we determine that 
A, the coefficient matrix is of rank L, and C, the augmented matrix, 
is of rank r. Then, 

Theorem 4 - The system AX = K has a solution if and only if L = r. 
Two possibilities exist: 

(1) L = r = n. In this case, D would be of the form 



and 



A,X 



1 

1 

9 ° 

1 o 
o 1 





A, = 



1 
1 

6 6 



) 




) 

















x x + o.x s + — + o.x n 
o.x, + x n + — + o.x„ 



O.X, + 0.X- + 
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from which we obtain X x = r x , , X n = r n . The system is consistent 

and has a unique solution. 

(2) L = r < n. An example of such a system is the following: 



for which 



C = 



and 



A = 





+ 2X a - 


x 3 


- 2X 4 


= -1 




+ X a + 


x 3 




x 4 


= 4 




- + 


2X 3 


+ 




= 5 




+ 3Xa - 


2X3 


- 3X 4 


- -3 


"l 


2 -1 


-2 


-1 






2 


1 1 


-1 


4 






1 


-1 2 


1 


5 






1 


3 -2 


-3 


"3. 







1 





1 





3 





1 


-1 


-1 


-2 

































"1 


2 


-1 


-2 




"l 





1 


-1 




2 


1 


1 


-1 







1 


-1 


-1 


1 


-1 


2 


1 
















_1 


3 


-2 


-3 













q 


2 


and 


L = 


2 


L = r < n 
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1 





1 







X' 




+ 


O.X, 


+ x 3 


+ 


O.X 4 = 3 







1 


-1 


-1 






= o.x t 


+ 


X3 


- x 3 




X 4 = -2 


















x 3 




+ 


o.x a 


o.x 3 


+ 


O.X 4 = 




















O.Xj 


+ 


o.x 3 


+ o.x 4 


+ 


O.X 4 = 



and 

X x = -X 3 + 3 

(1) 

X a = +X 3 + X 4 -2 

We have found X x and X, in terms of the remaining unknowns which 
can be considered parameters. 

The system is consistent and it is said to have a two-fold infinity 
of solutions; all solutions Xj , X a , X 3 , X 4 are formed by assigning 
arbitrary values to X 3 and X 4 , and determining X x and X, from (1), 
known as the general solution. 



On the other hand, if L < r s n, the system is inconsistent; that is, 
it has no solution. An example of such a system is the following: 

- Xj, + X 8 + X 3 = 1 

-5X 1 + 5X S + X 3 = 1 

X a - ^ + X 3 = 2 

for which 



-1 1 


1 


1 




1 


-1 





1/6 


-5 5 


1 


1 










1 


11/6 


J - 1 


1 


2 













-2/3 



and 





-1 


1 


1 




1 


-1 





A = 


-5 


5 


1 










1 




1 


-1 


1 
_ 


■ 
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By Theorem 3, r = 3, L = 2 /. L < r 





"l 


-1 


o" 






x x 


- X, 


+ 


o.x 3 


- 1/6 




^X = 








1 




x. 


- O.Xj 


+ o.x, 


+ 


x 3 


- 11/6 


(2) 









q 




. X 3. 


o.x x 


+ o.x, 


+ 


0.X3 


= -2/3 





But since the last equation in the set (2) cannot be satisfied by 
any (X x , X,, X 3 ) , we conclude that the system is inconsistent. 

Note that if the n x n system is homogeneous, it is consistent, and 
its solution depends on the value of L, the rank of the coefficient 
matrix. If L = n, the system has a unique solution; namely, the 
trivial solution. If L < n, the system has an infinite number of 
solutions . 

The method of reducing the augmented matrix C to its equivalent matrix 
D is carried out by operations on C registered in the same storage lo- 
cations where the original matrix C is placed at the beginning of the 
program. 

The total number of divisions and multiplications carried out for 
the reduction of a system containing a unique solution is given by 
the formula 

[ n 3 + (2L + 1) n a + 2n ] 

where 

n = number of equations 

L = total number of constant vectors 

When the reduction is completed, the matrix D is in storage, and an 
analysis of its elements will determine the correct solution. 
Obviously, the problem of computed zero values mentioned in the 
Program Description arises in this part of the problem. This problem 
and its solution is the subject of the next section of this program 
write-up. 

The time for solution of any system is approximated by 



T « a [V 3 +(2L + 1) n 8 + n "jp 
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LS= LEFT SIDE OF 

CHECK EQUATION 
RS=RIGTH SIDE CONSTANT 
OF CHECK EQUATION 
L=ROW RANK OF COEFF. 

MATR I X 
R=ROW RANK OF AUGMEN. 
MATR I X 
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where is the multiplication factor of the FORTRAN multiplication 
operation, and 1 s a £ 3 is a factor which takes care of additions, 
bookkeeping, and special routines if an infinite number of solutions 
exist. 

The procedure taken to produce the answers is shown schematically by 
the flow chart on page 10. This flow chart is a functional flow 
chart and, as such, does not show all the steps taken by the program, 
A detailed flow chart of the program is shown on pages 51 through 61. 



Determination of' the Zero Value 

As mentioned before in the Program Description and suggested by the 
theory explained in the Method of Computation, when checking for in- 
consistency the elements of the equivalent augmented matrix D must 
be compared against a calculated zero value. Also, it is desirable 
to estimate how accurate the calculated unknowns satisfy the system 
under solution. 

Unfortunately, the computer for which this program was set up consists 
of a 20K memory and uses a floating point arithmetic which carries 
eight significant digits only. This type of precision very s eld on re- 
sults in computed zero value and exact calculation of the unknowns in 
the system. 

For example, the matrix 
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31 


62 


33 


99 


27 


36 
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when reduced by computer operations showed the elements enclosed in 
dashed lines as 

- . 30000000E-06 . 40000000E-06 

.20000001E+01 .30000000E-06 
.60000000E-06 .29999997E+01 

. 14000000E-05 - .24400000E-05 

In order to take care of the above restriction, it was necessary to 
introduce a lower bound for the value of computed zero; this value 
does not only reveal possible zero elements , but it is also a measure 
of the accuracy carried through computations. This lower bound is 
determined by the relative size of the elements present in the aug- 
mented matrix, and by the way floating point arithmetic is handled in 
the computer itself. 

We have mentioned before that the reduction of the system's augmented 
matrix to its equivalent is carried out through a process of divisions, 
multiplications, and subtractions. In the course of the calculations, 
some elements become small with respect to others; and, in general, 
these elements are likely to contain larger percentage of errors; con- 
sequently, if division by any one of these elements takes place, the 
error is propagated to all elements in the matrix, at a fast rate. 
For example: 

Let X - .X 1 X a X3X 4 X 5 X 6 X 7 X B E+NM and Y = . Y 1 Y a Y 3 Y 4 Y s Y s Y 7 Y 8 E+NL 

be two elements appearing during the calculations . Let us assume 
that NL = NM-1, X-Y < 10, and Y 8 > 5. If the element Z = X-Y is 
calculated, the following happens 

.X 1 X 3 X 3 X 4 X 5 X 6 X 7 X S E+(NM) 
- .0 Y 1 Y 3 Y 3 Y 4 Y 5 Y a Y 7 E+(NL + 1 = NM) 

.0 Z 3 Z 3 Z 4 Z s Z 8 Z 7 Z a E+NM 

since X 8 is considered inaccurate, and Y 3 > 5 has been neglected, the 
Z 8 digit is also inaccurate. Furthermore, since Z 1 = 0, the number Z 
will be taken as Z = .Z s Z 3 Z 4 Z B Z e Z 7 Z 8 0E+(NL = MN-1), which carries the 
error in the seventh digit; division by this number will propagate this 
error . 
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The program, however, is set up so that when the operation division 
is to be carried out, it will be carried out in the row which con- 
tains the largest element (largest value indicates the largest number 
in absolute value); in this way, the "round-off" errors can be made 
as small as possible. 

Even though the above precaution has been taken, the following 
possibility exists. Suppose L = .L 1 L a L 3 L 4 L B L 6 L 7 L 8 E+LM is the 
largest element in the augmented matrix; depending on its position 
in the matrix, this element may be converted into another whose 
exponent may reach the value of LM+1; let us call this element 
L« = .L^L^L'L'I^L.L^ E+(LM+1) . 

If, during the calculations, A and B are two elements such that 
Z = L' - AB- = L' -C = 0, the following may happen 

.LjLjLs'LjL^LjLj E+(LM+1) 
- .C x C B C a C 4 C B C s C T C 8 E+(LM+1) 

.0 Z B E+(LM+1) 

the digit Z 8 appears in the answer, since L^ and C 8 are considered 
inaccurate according to FORTRAN manuals . 

Knowing that this is the worse case that can appear during calculation, 
the value Z = .Z a 0000000 E+(LM+l-7) = .Z 8 0000000 E+(LM-6) was assumed 
to be a good "lower bound" for the value of zero. Z a could be anywhere 
between one and nine, but actual calculations showed that Z 8 = 1 was 
accurate enough. 

This value Z is considered to be the zero of the system, and any element 
whose absolute value becomes smaller than or equal to Z is considered a 
zero element by the program. 

For any (n x n) system, Z = . 1000000E(EX) , where EX is the number 
resulting from the subtraction of six from the exponent of the 
largest element in the augmented matrix. Thus, if 350 is the lar- 
gest element, Z = .1000000E-03. 

This value of Z must be entered into the program from the typewriter, 

after the message 

ENTER TOLERABLE Z IN E14.8 FORMAT SPEC. 



is typed . 
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Z is a greatest lower bound and should not be selected any smaller. 
On the contrary, if it seems to be very small for the amount of 
operations carried during a given calculation, a message is typed 
indicating that Z should be increased in value. 

This zero value also implies that the relative size of the elements 
forming the augmented matrix, shall not differ by more than a factor 
of 10 e or shall not be any closer than a factor of 10" s . Otherwise, 
the accuracy called for goes beyond the assumed zero value, which is 
a contradiction. Such systems shall be considered as " ill-conditioned " 
systems and are out of the scope of this program. 

One may argue that the system 

10 7 Xj + 10 7 Xg = 10- 3 

2 x 10 7 X x + 3 x 10 7 X a = 10- 3 

is very consistent, and its unknowns are X x = -7 x 10" 10 and 
X 3 = 8 x 10- 10 . However, if analyzed by this program, Z = 
•10000000E+02, and the answers could be interpreted as zeroes. 

But this is not quite an ill-conditioned system, since it could be 
written as : 

X x ' + X s ' = 10" 3 XI = 10 7 X 1 , X^ = 10 7 ^ 

2XJ + 3x a' " 10-3 

for which Z = .1000000E-05 and X x * --7 x 10" 3 , X^ = 8 x 10 -3 , which 
implies X x = -7 x 10 _1 ° , and Xg = 8 x 10 -10 

The system 

10 s X 1 + 10 s X a = 10 " s 
.34 X x + 10 s X 8 = 10 7 

for which Z = .10000000E+02 is a true ill-conditioned system, and 
could not be handled by this program, not only because of the zero 
value chosen, but because the numbers are so far apart that single 
precision arithmetic could lead to very inaccurate answers. 

Such problems could be solved if the program was set up for multiple 
precision arithmetic. 
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Input Card Format 

This program is set up to solve systems of N equations with N 
unknowns only. 

If the system contains L equations and N unknowns and L < N, a 
number of N-L equations should be added after the L equation, be- 
fore attempting to use this program with such a system. These 
extra equations must consist of zero coefficients and zero constant 
equalities. For example, the system 



shall be transformed i 







+ 


a ia 


X, 


+ 


a !3 


x 3 


- 




X x 


+ 


a aa 


X a 


+ 


a 3 3 


x 3 


= 


Lnto 


















a n 


Xi 


+ 


a i 8 


x a 


+ 


*13 


x 3 


= 




Xx 


+ 




x a 


+ 


a 3 3 


x 3 


- 


0. 


XI 


+ 


0. 


Xx 


+ 


0. 


x 3 


- 



For any N x N system (smaller than or equal to 17 x 17) one can obtain 
solutions for as many as 18-N constant vectors, which will carry the 
names "Vector N+l, Vector N+2, , Vector 18." 

The first card in the input data must contain two numbers, N and M. 
N is the number of unknowns (or the number of equations) forming the 
system, and M is the number of unknowns plus the number of constant 
vectors. These two numbers must be punched in the FORMAT (13, 13) 
specification form; that is, the units position of number N must be 
column 3, and the units position of number M must be column 6. 

The rest of the cards must contain the elements of a rectangular 
matrix formed by the system coefficient matrix and all the constant 
vectors added as columns to the right of it. These elements must be 
punched row wise starting with the element in the first column and 
ending with the element in the last column. Each card must contain 
six elements, and they must be punched in the E12.6 FORMAT specifica- 
tion form. Depending on the number of elements per row, M to be 
exact, a minimum of one card, or a maximum of three cards per row will 
exist. If M is not a multiple of six, zero values must be introduced 
in order to complete six numbers in the last card of a row. A number 
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of N rows must be present in the input. 1 Thus, depending on the 
value of M, the following table gives the total number of cards 
present in the input: 



Value of M 


Number of Cards 


M s. 6 


N + 1 


6 < M < 12 


2N + 1 


12 < M £ 13 


3N + 1 



For those unfamiliar with the FORMAT specifications, the specification 
E12.6 means the following: The number under this specification must 
take a total of 12 spaces. Space one is for the sign of the number, 
space two for the decimal point, spaces three through eight for the 
significant digits of the number and zeros if required, and spaces 
nine through twelve for a minus sign if the exponent is negative or 
a plus or blank if it is positive, and two spaces for the exponent. 
The following examples show how the numbers on the left must be punched 
on the card: 



Number 



Punched 



35 +.350000E+02 or b.350000E+02 

785.25 b.785250E+03 
-.00032 -.32O000E-03 

b indicates a blank space 



Output Format 

The output from this program is printed by means of the console 
typewriter, and consists of four types: 



Note: If the first column of the coefficient matrix consists of 
zero elements only, the system should be rearranged so that this 
column appears in a different position. 
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(1) The message "RELOAD DATA INCR. Z BY FACTOR OF 10." This 
message is typed to indicate that the computed unknowns do not satisfy 
the system within the accuracy of the tolerable Z value. 

After this message is typed, the operator should reload 
the data and push the start button. Once the program requests the 
entering of the tolerable Z value, the operator shall enter a new 
Z value according to the corresponding instructions found under 
Operating Instructions. 

This new Z value is determined from the equation 

Z - 10 x Z x 

where Z x is the tolerable Z value used previously. 

(2) If the system is consistent, and the answers computed 
satisfy the accuracy assigned by Z, the message 

SOLUTION FOR VECTOR N+l 

is typed first, to indicate that the column of numbers immediately 
below it are the answers corresponding to the first constant vector. 
The answers are next typed one per line, in the E14.8 FORMAT specifi- 
cation, and they correspond to X x , X,...., X fl . 

If more than one constant vector has been included in the 
input data, the program continues to type SOLUTION FOR VECTOR N+2, 
and will continue as above until the solution for the last vector 
(vector M) has been typed. 

(3) If together with any one of the constant vectors, let us 
say vector N+l, the system is inconsistent, the message 

SOLUTION FOR VECTOR N+l 
SYSTEM INCONSISTENT 

is typed . 

(4) If together with any one of the constant vectors, let us 
say vector M, the system has an infinite number of solutions, the 

message 

SOLUTION FOR VECTOR M 
INF SOLS SW3 ON FOR GEN S0L-SW3 OFF TO CONT. 

is typed, and the program halts. 
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This message tells the operator that the system has an 
infinite number of solutions, and it gives him the choice of obtaining 
a general solution or continuing with the solution for the next vec- 
tor, or for the next system, whichever happens to be next. 

If the operator does not want to obtain a general solution, 
he should set SW3 OFF and PUSH START. 

If the operator wishes to obtain a general solution, he 
should set SW3 ON and PUSH START. 

The output will then consist of (N-L+l) columns, each 
containing L numbers in the E14.8 FORMAT specification; L is the 
row rank of the coefficient matrix. The first (N-L) columns are 
headed by the heading (1 < i <. n) , and the last column appears 
without any heading, but reasonably spaced from the column before. 

The Xi's indicate those variables that can be considered 
as parameters, and the numbers immediately below them indicate their 
coefficients in the general solution. The last column consists of 
the constants appearing in the general solution. 

How to form a general solution as (1) in the Method of 
Computation section is better explained through an example. 



The system 



3X X 


+ 


3X, 


+ 


4X 3 


+ 


8X 4 


+ 


5X 5 


+ 


15X 6 


+ 


9X 7 


= 12 


2Xj 


+ 


2X S 


+ 


3X 3 


+ 


6X 4 


+ 


7X S 


+ 


21X 6 


+ 


6X 7 


= 8 


Xi 


+ 


x 8 


+ 


5X 3 


+ 


10X 4 


+ 


8X 5 


+ 


24X 6 


+ 


3X 7 


= 4 


4X X 


+ 


4X 3 


+ 


10X 3 


+ 


20X 4 


+ 


11X S 


+ 


33X a 


+ 


12X 7 


- 16 


6X 1 


+ 


6X 8 


+ 


IIX3 


+ 


22X 4 


+ 


13X 5 


+ 


39X 6 


+ 1.8X 7 


= 24 


5^ 


+ 


5X 3 


+ 


I2X3 


+ 


24X 4 


+ 


12X 5 


+ 


36X 6 


+ 


15X 7 


= 20 


9X, 


+ 


9X 3 


+ 


31X 3 


+ 


62X 4 


+ 


33X 5 


+ 


99X 6 


+ 


27X 7 


= 36 



for which Z = .1000G000E-04 produces the following answers 
RELOAD DATA INCR. Z BY FACTOR OF 10 
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Once the data is reloaded and Z is made .0000000E-03, 
the following answer is obtained: 

SOLUTION FOR VECTOR 8 
INF SOLS SW3 ON FOR GEN S0L-SW3 OFF TO CONT. 

If SW3 is ON, the following is obtained: 

X a 

- . 10000000E+01 
.00000000E-99 
.000O000OE-99 
.OOOOOOOOE-99 



X 4 

.30000000E-06 
-.20000001E+01 
-.60000000E-06 
-.14000000E-05 



.40000000E-06 
.30000000E-06 
.29999997E+01 
.24400000E-05 

.40000000E+01 
.OOOOOOOOE-99 
.OOOOOOOOE-99 
.00000000E-99 



Since Z was originally estimated as .10000000E-04, all of those 
numbers whose exponent is smaller than E-04 are considered zero. 

From above, form a table consisting of four columns 
(N-L+l in the general case), and four rows (L in the general 
case). The first column in this table must consist of the four 
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numbers appearing under , and X a is its heading; the second column 
consists of the numbers under X 4 , and X 4 is its heading; the third 
column consists of the four numbers under Xg , and X e is its heading; 
the fourth column consists of the four numbers which appear in the 
column without heading. Thus so far, we have 



x 8 


x 4 


x s 




-1 








4 





-2 














-3 


















Now, on the left side, assign names to the rows, by means of the 
remaining unknowns Xj , X 3 , X B , X 7 , and in the same order. Thus, 
at the end, the table looks like 





x 8 


x 4 


Xe 






-1 








4 


x 3 





-2 








x 5 








-3 





X, 















and the general solution is obtained as 

X x = -l.Xg + O.X 4 + O.X 6 +4 

X 3 = O.X a - 2.X 4 + O.X 6 + 

X e - C.X 2 + O.X 4 - 3.X 6 + 

X 7 = O.X 8 + O.X 4 + O.X 6 + 
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or 

X, = -X a + 4 

X 3 = -2X 4 X, , X 4 , X, take on any value 

X 5 = -3Xg 

X 7 = 
Operating Instructions 
I . A. Initial Console Setting 





Program 


Stop 


Parity 




X 


I/O 




X 


Overflow 


X 





B. Sense Switch Settings 



sense switch 
1 


OFF 


Mot used 


Sense Switch 
2 


OFF 


Not used 


Sense Switch 
3 


ON 
OFF 


If inf. Sols, exist, general solution is typed 
General solution is not typed 


Sense Switch 
4 


ON 
OFF 


To correct error in typing tolerable zero value 
Tolerable zero value entered correctly 



See section V for further comment on switches #3 and #4. 
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II . Input-Output 

Card Reader 



No. of Cards 


Description 


rn-i 

joZ. 


Program Deck 


1 


N - number of equations in first set of 
N&M Card - equations 

M - number of equations plus number of con- 
stant vectors - first set of equations 


2-51 


Coefficient cards for first set of equations 


1 


„.„ , N - number of equations in second set of 
N&M Card - ^ 
equations 

M - number of equations plus number of con- 
stant vectors - second set of equations 


2-51 


Coefficient cards for second set of equations 




ETC. 



Typewriter Output 



Control 


Description 


Margins 


Set left margin as desired 


Tab Stops 


None 


Forms 


Standard paper 



Page 23 



III . Normal Loading Procedure 

(1) Clear Storage 

(2) Depress RESET 

(3) Place the PROGRAM DECK 1 in the reader hopper and depress 
LOAD button. When the READER NO FEED light comes on, depress the 
READER START so that the last two program cards are read in. At 
this time, the message LOAD DATA is typed. 

(4) Place DATA cards in the reader hopper and depress COMPUTER 
START; when the READER NO FEED light comes on, depress the READER 
START so that the data cards are read into memory. 

If the data consists of one system of equations only, the 
READER NO FEED light will come on, at which time the READER START 
must be depressed to read the last two data cards . If the data con- 
sists of two or more systems, the READER NO FEED light will come on 
when the computer is attempting to read the data corresponding to the 
last system; again, the READER START must be depressed to read the 
last two data cards for this system. 

IV. Special Loading Instruction 

After the program deck has been read into memory and it has been 
started, it can be stopped at any time by depressing COMPUTER INSTANT 
STOP. 

The program can then be initialized and started by depressing the 
RESET button and inserting the instruction 4908000 by means of the 
console typewriter. 

V . Special Instructions and Remarks 

Tolerable Zero Value : When the program begins, it first reads the 
data for the first system, and then the message 

ENTER TOLERABLE Z IN E14.8 FORMAT SPEC. 

is typed. The operator must then type the Z value that applies for 
the system under analysis, and in accordance with the E14.8 FORMAT 
specification. That is, the Z value must be entered as +.10000000E(EX), 
where EX is the exponent determined as explained before in the "Determina- 
tion of Zero Value" section. 



1 Recall the definition appearing on page 1 where the Program Deck is 
defined as Deck #2 without the hash total card; or, it can be obtained 
by compiling the S ource Program (Deck #1 without hash total card) by 
means of a FORTRAN 1 with FORMAT compiler containing the Absolute value 
subroutine . 
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If the Z value is typed correctly, set SW4 OFF and depress RELEASE 
and COMPUTER START. The program will then continue to solve the 
system under analysis. 

If the Z value is not typed correctly, set SW4 ON and depress RE- 
LEASE and COMPUTER START. The program will immediately return 
control to the typewriter so that the Z value may be entered again. 
If the Z value is re-entered correctly, follow the procedure for a 
correct entry. 

If the Z value is typed incorrectly, and SW4 is OFF, depressing the 
RELEASE and START buttons will produce the message ERROR F7 and the 
program will continue. Continuation of the program Under these cir- 
cumstances is useless; consequently, the programmer must press the 
COMPUTER INSTANT STOP. 

After the program is stopped, the programmer must place the non-pro- 
cessed system data cards back in the reader hopper. He does this in 
the following way: 

(a) Remove data cards, if any, from the reader hopper. 

(b) Depress NON-PROCESS RUN-OUT key on the card reader. 

(c) Remove cards, if any, from the error select stacker. 

(d) Place these cards in front of the cards removed from 
the hopper. 

(e) Remove non-processed system data cards from the non- 
select stacker. 

(f) Place these cards in front of the cards removed from the 
error-select stacker and replace deck in reader hopper. 

With the non-processed system data cards ready for processing, the 
programmer must depress the COMPUTER RESET and initialize the pro- 
gram as explained in section IV. When the READER NO FEED light comes 
on, the READER START is depressed and the program will continue as 
explained in previous sections. 

Accuracy Warning : As explained under Output Format section of this 
write-up, when the computed unknowns do not satisfy the system under 
analysis within the accuracy of the tolerable Z value, the message 



RELOAD DATA INCR. Z BY FACTOR OF 10 
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is typed; after this the program halts waiting for the programmer's 
action. 

If the programmer does not want to continue with this system, but 
he wants to continue with the next system for solution, he does so 
by simply depressing the COMPUTER START. 

If the programmer wants to continue with the solution of the system 
for which the warning message was typed, he does so by reloading the 
system's data cards as explained in this section under the title 
Tolerable Zero Value. 

Once the data cards are in the hopper, pushing the COMPUTER START 
and the READER START will initialize the program. When the request 
for a new zero value is typed, the programmer follows the instruc- 
tions insinuated in the"0utput Format" section of this write-up. 

General Solution : When the system under analysis has an infinite 
number of solutions, the message 

INF SOLS SW3 ON FOR GEN S0L-SW3 TO CONT. 

is typed; after this the program halts waiting for the programmer's 
action. 

If the programmer does not want to obtain a general solution, he sets 
SW3 OFF and depresses START. The program will continue with the solu- 
tion for the next constant vector or the next system, if any. If the 
program tries to continue with the next system, but there is no data 
cards for a next system, the program stops on a READER NO FEED warning 

If the programmer wants to obtain a general solution, he does so by 
setting SW3 ON and depressing START. 

VI . Programmed Stops and Required Action 

All programmed stops are accompanied by typewriter messages which are 
self-explanatory and which indicate the required action. 

As the program is written in FORTRAN, all the error stops and 
messages of this system apply. 
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Sample Problems 

The following seven sample problems were used one after the other 
and they illustrate the program reliability: 



irroblem No. 1: 





+ 


2Xa 


- x 3 


- 2X 4 = 


-1 


2X, 


+ 


x 3 


+ X 3 


- x 4 = 


4 


Xi 




x 8 


+ 2X 3 


+ x 4 = 


5 


X i 


+ 


3X 8 


- 2X3 


- 3X 4 - 


-3 



Z = .10000000E-G5 



SOLUTION FOR VECTOR 5 
INF SOLS SW3 ON FOR GEN S0L-SW3 OFF TO CONT. 



.10000000E+01 
.lOOOOOOOE+01 

X 4 

.00000000E-99 
.10000000E+01 

.30000000E+01 
.20000000E+01 





x 3 


x 4 




X, 


X! 


-1 





3 


x 8 


x a 


1 


1 


-2 


x 3 



on any rea 
value 



Problem No. 2: 



3X X + 8X3 + 6X3 + 10X 4 + 42X 5 = 

2X3 + X 4 + 5X 5 = 

.OIX3 + 4X 4 + 6X B = 

2X, + 4X + 7X, + 9X„ = 



Z = .10000000E-0 
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Add one more equation, as explained in the Input Card Format 
section, so that the system becomes 

3Xj + 8X 3 + 6X3 + 10X 4 + 42X 5 - 

2X 3 + + X 4 + 5X 6 = 

.OIX3 + 4X 4 + 6X B = Z = .10000000E-04 

2Xj + 4X 3 +7X3 + 9X 6 = 

O.X a + O.Xg + O.X3 + O.X 4 + O.X E = 

SOLUTION FOR VECTOR 6 
INF SOLS SW3 ON FOR GEN S0L-SW3 OFF TO CONT . 



.87446120E+01 
.17472361E+01 
.22111665E+01 
.15055279E+O1 

.00000000E-99 
.00000000E-99 
.OOOOOO00E-99 
.00000000E-99 



roblem No. 3: 





x s 




Xi 


-8.74461 





x 3 


-1.74724 





x 3 


2.21117 





x 4 


-1.50553 






X x - -8.74461 X 5 

Xa - -1.74724 X s 

X 3 - 2.21117 X 5 

X 4 = -1.50553 X fi 

X 5 = any real 
constant 



X x 


+ 




+ 


3X 3 


- 4X 4 


+ 5X 5 - 


6X 8 


+ 


7X 7 


- 8X 8 - 


-28 


-14 


-84 


60 






7X 3 








- 2X 5 




+ 


X 7 


+ x 8 = 


19 


9.5 


57 


30 




+ 


3X 2 




4X3 


- 5X 4 




6X 8 








12 


6 


36 


90 


X l 


+ 


x a 


+ 


x 3 


+ x 4 


+ X 5 + 


x 6 


+ 


X 7 


+ X 8 = 


36 


18 


108 


45 


10X l 


+ 


3X a 




4X 3 




+ X B + 


2X a 




9X 7 




-42 


-21 


-126 


135 


3X X 




3X, 




2X 3 


+ 2X* 


+ 


x 8 








5 


2.5 


15 


67.5 


-8X l 








9X 3 




+ 7X 5 




+ 


6X 7 


- 3X 8 = 


18 


9 


54 


99 




















3X 7 


- 5X e = 


-19 


-9.5 


-57 


































































Z = .10000000E-03 



SOLUTION FOR VECTOR 9 
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.99999970E-00 
.1999999SE+01 
.29999999E+01 
.40000002E+01 
.49999990E+01 
.60000010E+01 
.69999996E+01 
.79999996E+01 



X, 



= 1 

= 2 

= 3 

= 4 

= 5 

= 6 

= 7 



SOLUTION FOR VECTOR 10 
X, = .5 Xg = 1.0 X3 = 1.5 X 4 = 2.0 
X 5 = 2.5 X a = 3 X 7 = 3.5 X B = 4 

SOLUTION FOR VECTOR 11 
Xj = 3 X 3 = 6 X 3 = 9 X 4 = 12 
X 5 =15 X 6 = 18 X 7 =21 X 8 = 24 



X, - 



SOLUTION FOR VECTOR 12 

.17373927E+02 X 4 = .23638003E+01 X 7 = .12568545E+02 

.46417256E+01 X 5 = .11300876E+02 X a = .75411273E+01 
-.11788599E+02 X R = .99859570E-00 



Problem No. 4: 

The system whose coefficient matrix is shown on page 29 will have 
an infinite number of solutions or will be inconsistent, depending 
on the values of the constant vectors. This is so because the 
determinant of the system is zero. 



This coefficient matrix, together -with the constant vectors shown 
to the right and Z = .10000000E-03, gives the following answer: 
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SOLUTION FOR VECTOR 14 

SYSTEM INCONSISTENT 
SOLUTION FOR VECTOR 15 

SYSTEM INCONSISTENT 
SOLUTION FOR VECTOR 16 

SYSTEM INCONSISTENT 
SOLUTION FOR VECTOR 17 

SYSTEM INCONSISTENT 
SOLUTION FOR VECTOR 18 

INF SOLS SW3 ON FOR GEN SOL-3W3 OFF TO CONT. 






35 


10 


40 


10 


30 


150 


115 


40 


41 


53 


55 


67 




141 


17 


20 


38 





- 35 





11 


15 


17 


13 


12 


61 


64 


63 


72 


75 


47 




48 


49 


52 


15 


- 70 


- 10 


-11 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 




11 


12 


13 


14 


- 20 


- 40 


-15 


- 1 





11 


15 


32 


23 


17 


11 


15 


16 


17 




18 


19 


20 


21 


- 80 


- 10 


-17 


- 2 


-11 





150 


149 


148 


147 


146 


145 


144 


143 




142 


141 


140 


139 


- 20 


- 30 


-13 


- 3 


-15 


-150 





- 2 


- 3 


- 4 


-5 


- 6 


- 7 


- 3 




- 9 


-10 


-11 


-12 


- 60 


-150 


-12 


- 4 


-32 


-149 


2 





21 


22 


23 


24 


25 


26 




27 


28 


29 


30 


-300 


-115 


-61 


- 5 


-23 


-148 


3 


-21 





31 


35 


39 


43 


47 




147 


148 


149 


150 


-230 


- 40 


-64 


- 6 


-17 


-147 


4 


-22 


-31 





160 


250 


349 


471 




582 


651 


781. 


811 


- 80 


- 41 


-63 


- 7 


-11 


-146 


5 


-23 


-35 


-160 





-30 


- 40 


-50 




-60 


-70 


-80 


-90 


- 82 


- 53 


-72 


- 8 


-15 


-145 


6 


-24 


-39 


-250 


30 





- 13 


-15 




-14 


-16 


-15 


-17 


-106 


- 55 


-75 


- 9 


-16 


-144 


7 


-25 


-43 


-349 


40 


13 





111 




121 


133 


144 


155 


-110 


- 67 


-47 


-10 


-17 


-143 


8 


-26 


-47 


-471 


50 


15 


-111 







_ 11 


9 


7 


5 


-134 
_ 



COEFFICIENT MATRIX CONSTANT VECTORS 



EXAMPLE 4 
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-.89912030 x X 13 


x 3 


= 


-1.1049142 x X 13 


x 3 


- 


57.970323 x X 13 


x 4 


= 


-9.417142 x X 13 


x s 


= 


.14746840 x X 13 


x 6 




1.7722926 x X 13 


x 7 




7.5299710 x X 13 


x 8 




-12.344710 x X 13 


x* 




1.176188 x X l3 






6.271787 X 13 






-3.5682707 X 13 


x 13 




-2.0166340 X 13 



X 13 could take on any value 
Problem No. 5 : 

The system whose coefficient matrix and constant vectors are shown 
on page 31 and for which Z = .10000000E-03 produces the following 
answers : 

RELOAD DATA INCR . Z BY FACTOR OF 10 

with Z = .10000000E-02 

SOLUTION FOR VECTOR 14 SOLUTION FOR VECTOR 15 

-.89017740E+01 - .77705800E+01 

-.48386103E+01 - .49968307E+01 

.41425661E+03 .37400470E+03 

5548039 1E+02 - .51928873E+02 

.14276311E-00 .33781753E-00 

.15499623E+02 .15138997E+02 

.56783579E+02 .49755229E+02 

-.99629510E+02 - .88164190E+02 

.10310057E+02 .88469110E+01 
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35 


10 


40 


10 


30 


150 


115 


40 


41 


53 


55 


67* 




141 


17 


20 


38 





- 35 





11 


15 


17 


13 


12 


61 


64 


63 


72 


75 


47 




48 


49 


52 


15 


-70 


- 10 


-11 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 




11 


12 


13 


14 


-20 


- 40 


-15 


- 1 





11 


15 


32 


23 


17 


11 


15 


16 


17 




18 


19 


20 


21 


-80 


- 10 


-17 


- 2 


-11 





150 


149 


148 


147 


146 


145 


144 


143 




142 


141 


140 


139 


-20 


- 30 


-13 


- 3 


-15 


-150 





- 2 


- 3 


- 4 


- 5 


- 6 


- 7 


- 8 




- 9 


-10 


-11 


-12 


-60 


-150 


-12 


- 4 


-32 


-149 


2 





21 


22 


23 


24 


25 


26 




27 


28 


29 


30 


-300 


-115 


-61 


- 5 


-23 


-148 


3 


-21 





31 


35 


39 


43 


47 




147 


148 


149 


150 


-230 


- 40 


-64 


- 6 


-17 


-147 


4 


-22 


-31 





160 


250 


349 


471 




582 


651 


781 


811 


- 80 


- 41 


-63 


- 7 


-11 


-146 


5 


-23 


-35 


-160 





-30 


- 40 


-50 




-60 


-70 


-80 


-90 


- 82 


- 53 


-72 


- 8 


-15 


-145 


6 


-24 


-39 


-250 


30 





- 13 


-15 




-14 


-16 


-15 


-17 


-106 


- 55 


-75 


- 9 


-16 


-144 


7 


-25 


-43 


-349 


40 


13 





111 




121 


133 


144 


155 


-110 


67 


-47 


-10 


-17 


-143 


8 


-26 


-47 


-471 


50 


15 


-111 







11 

- 


9 


7 


5 


-134 



COEFFICIENT MATRIX CONSTANT VECTORS 



EXAMPLE 5 
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•47450099E+02 
-.12369882E+02 
-.26474642E+02 

.91686745E+02 



.39042339E+02 
-.89736869E+01 
-.23229999E+02 

.81447152E+01 



SOLUTION FOR VECTOR 16 SOLUTION FOR VECTOR 17 SOLUTION FOR VECTOR 18 



-.70374340E+01 
-. 44552 765E+01 
.33750718E+03 
-.46798547E+02 
.29713861E-00 
.13666631E+02 
. 45032 153E+02 
-.79983800E+02 
.79584935E+01 
.34742553E+02 
-.74102383E+01 
-.20863395E+02 
.74716143E+01 



-.68729740E+01 
-.50068592E+01 
.34873410E+03 
-.50506716E+02 
.45286507E-00 
.12175005E+02 
. 4773208 7E+02 
-.82814740E+02 
.84018230E+01 
.38005589E+02 
-.95950557E+01 
-.2 140448 5E+02 
.79367011E+01 



.96000000E-05 
-.24577692E+01 
. 128948 79E+03 
-.20947452E+02 
.32802864E-00 
•39422583E+01 
. 16749622E+02 
-.27457473E+02 
.261630O2E+01 
. 13950905E+02 
-. 7937243 7E+01 
-.44857705E+01 
.22243907E+01 



Note that the value of X 1 in the solution for vector 18 is smaller 
than Z - .10000000E-03 and consequently X, =0; this is the true 
value of Xj in the system. 

Problem No. 6 ; 

The system shown on page 33 whose determinant is clearly zero and 
for which Z = .10000000E-03, produces the following answers: 

RELOAD DATA INCR. Z BY FACTOR OF 10 
with Z = .10000000E-02, we obtain 

X t - -.48874995 X^ + 2.000 
X, - 2.397§300 X 17 
X 3 = -28.744234 X 17 
X 4 - 9.7720987 X 17 
X 5 = -.61967767 X 17 



X. =-1.7527315 X l7 



n, 
U 


35 


IU 




1U 


jU 






40 


- 35 





11 


15 


17 


13 


12 


61 


64 


- 10 


-11 





1 


2 


3 


4 


5 


6 


- 40 


-15 


- 1 





11 


15 


32 


23 


17 


- 10 


-17 


- 2 


-11 





150 


149 


148 


. 147 


- 30 


-13 


- 3 


-15 


-150 





- 2 


3 


- 4 


-150 


-12 


- 4 


-32 


-149 


2 





21 


22 


-115 


-61 


- 5 


-23 


-148 


3 


-21 





31 


- 40 


-64 


- 6 


-17 


-147 


4 


-22 


- 31 





- 41 


-63 


- 7 


-11 


-146 


5 


-23 


- 35 


-160 


- 53 


-72 


- 8 


-15 


-145 


6 


-24 


- 39 


-250 


- 55 


-75 


- 9 


-16 


-144 


7 


-25 


- 43 


-349 


- 67 


-47 


-10 


-17 


-143 


8 


-26 


- 47 


-471 


-141 


-48 


-11 


-18 


-142 


9 


-27 


-147 


-582 


- 17 


-49 


-12 


-19 


-141 


10 


-28 


-148 


-651 


- 20 


-52 


-13 


-20 


-140 


11 


-29 


-149 


-781 


- 38 


-15 


-14 


-21 


-139 


12 


-30 


-150 


-811 



COEFFICIENT MATRIX 
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41 


53 


55 


67 


141 


17 


20 


38" 




r ~i 




63 


72 


75 


47 


48 


49 


52 


15 




- 70 


7 


8 


9 


10 


11 


12 


13 


14 




- 20 


11 


15 


16 


17 


18 


19 


20 


21 




- 80 


146 


145 


144 


143 


142 


141 


140 


139 




- 20 


- 5 


- 6 


- 7 


- 8 


- 9 


-10 


-11 


-12 




- 60 


23 


24 


25 


26 


27 


28 


29 


30 




-300 


35 


39 


43 


47 


147 


148 


149 


150 




-230 


160 


250 


349 


471 


582 


651 


781 


811 




- 80 





-30 


- 40 


-50 


-60 


-70 


-80 


-90 




- 82 


30 





- 13 


-15 


-14 


-16 


-15 


-17 




-106 


40 


13 





111 


121 


133 


144 


155 




-110 


50 


15 


-111 





11 


9 


7 


5 




-134 


60 


14 


-121 


-11 





15 


20 


30 




-282 


70 


16 


-133 


- 9 


-15 





1 


20 




- 34 


80 


15 


-144 


- 7 


-20 


- 1 





10 




- 40 


90 


17 


-155 


- 5 


-30 


-20 


-10 







- 76 



C.V. 18 



EXAMPLE 6 
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x 7 






x e 




. uzizyyyo x 17 


x e 




-.ouoo/iy x 17 


X 10 




1 Q "7 £ "7 n Y 


Xn 




3.5436282 X l7 


X 18 




1.9503036 X 17 


X 13 




.43309011 X l7 






-1.8223183 X l7 


x 1B 




10.547266 X l7 






-10.317607 X l7 


Xj 7 


could take on any value 



Problem No. 7 : 

The system shown on page 35 , for which Z = .10000000E-03, gives 
RELOAD DATA INCR. Z BY FACTOR OF 10 

Z = .10000000E-02 

RELOAD DATA INCR. Z BY FACTOR OF 10 

Z = .10000000E-01 

SOLUTION FOR VECTOR 18 

Xi = .91800000E-04 Xg = - .11702277E+01 

Xg = .11183719E+02 X 10 = - .87497130E+01 

X 3 - -.13360065E+O3 X 1X = . 16529699E+02 

X 4 - .45581915E+02 X 13 = .90972961E+O1 

X 5 = -.28904901E+O1 X 13 = .20201331E+01 

X s = -.81754060E+01 X 14 - - .85001620E+01 

X 7 - -.38309755E+01 X 1S - .49197987E+02 

X 8 = .99337500E-01 X 1B = - .48126945E+02 

X 17 = .46646303E+01 

Note that X 1 < .10000000E-03 , and consequently X x = 0; 
this is the true value of X x in the system 
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35 


10 


40 


10 


30 


150 


115 


40 


41 


53 


55 


67 


141 


17 


20 


38* 




0" 


- 35 





11 


15 


17 


13 


12 


61 


64 


63 


72 


75 


47 


48 


49 


52 


15 




- 70 


- 10 


-11 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 




- 20 


- 40 


-15 


- 1 





11 


15 


32 


23 


17 


11 


15 


16 


17 


18 


19 


20 


21 




- 80 


- 10 


-17 


- 2 


-11 





150 


149 


148 


147 


146 


145 


144 


143 


142 


141 


140 


139 




- 20 


- 30 


-13 


- 3 


-15 


-150 





- 2 


- 3 


4 


-5 


- 6 


- 7 


- 8 


- 9 


-10 


-U 


-12 




- 60 


-150 


-12 


- 4 


-32 


-149 


2 





21 


22 


23 


24 


25 


26 


27 


28 


29 


30 




-300 


-115 


-61 


- 5 


-23 


-148 


3 


-21 





31 


35 


39 


43 


47 


147 


148 


149 


150 




-230 


- 40 


-64 


- 6 


-17 


-147 


4 


-22 


- 31 





160 


250 


349 


471 


582 


651 


781 


811 




- 80 


- 41 


-63 


- 7 


-11 


-146 


5 


-23 


- 35 


-160 





-30 


- 40 


-50 


-60 


-70 


-80 


-90 




- 82 


- 53 


-72 


- 8 


-15 


-145 


6 


-24 


- 39 


-250 


30 





- 13 


-15 


-14 


-16 


-15 


-17 




-106 


- 55 


-75 


- 9 


-16 


-144 


7 


-25 


- 43 


-349 


40 


13 





111 


121 


133 


144 


155 




-110 


- 67 


-47 


-10 


-17 


-143 


8 


-26 


- 47 


-471 


50 


15 


-111 





11 


9 


7 


5 




-134 


-141 


-48 


-11 


-18 


-142 


9 


-27 


-147 


-582 


60 


14 


-121 


-11 





15 


20 


30 




-282 


- 17 


-49 


-12 


-19 


-141 


10 


-28 


-148 


-651 


70 


16 


-133 


- 9 


-15 





1 


20 




- 34 


- 20 


-52 


-13 


-20 


-140 


11 


-29 


-149 


-781 


80 


15 


-144 


- 7 


-20 


- 1 





10 




- 40 


38 


-15 


-14 


-21 


-139 


12 


-30 


-150 


-811 


90 


17 


-155 


- 5 


-30 


-20 


-10 







- 76 
— 



COEFFICIENT MATRIX C.V. 18 
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SAMPLE PROBLEMS- INPUT 



Problem No. 1 



4 5 

. 100000E+01 .200000E+01-. T00000E+01 -.200000E+01 -. 100000E+01 .OOOOOOE-99 

.200000E+01 .100000E+01 . 1 00OO0E+0 1 -. 1 00000E+01 .400000E+01 .OOOOOOE-99 

. 100000E+01-.100000E+01 .200000E+01 . 100000E+01 .500000E+01 .OOOOOOE-99 

. 1 00000E+0 1 . 300000E+0 1 -. 200000E+0 1 -. 300000E+0 1 -. 300000E+0 1 . OOOOOOE-99 



Problem No. 2 



5 6 

.300000E+01 
.000000E-99 
.000000E-99 
.200000E+01 
.000000E-99 



.800000E+01 

.200000E+01 

.000000E-99 

.400000E+01 

.000000E-99 



.600000E+01 
.000000E-99 

.100000E-01 

.700000E+01 
.000000E-99 



. 1 00000E+02 
. 100000E+01 
.400000E+01 

.000000E-99 
.ooooooE-99 



.420000E+02 

.500000E+01 

.600000E+01 

.900000E+01 
.000000E-99 



.000000E-99 
.000000E-99 
.000000E-99 
.000000E-99 
.000000E-99 



Problem No. 3 



8 12 

, 100000E+0.1 

.700000E+01- 
. 000000E-99 

. 100000E+01 
.200000E+01 

.000000E-99 

. 100000E+01 
.100000E+01 
. 1 00000E+02 

.900000E+01 

.300000E+01- 

. 000000E-99 

.800000E+01 

,6oooooe+oi- 

. OOOOOOE-99 

.300000E+01- 



.200000E+01 

.800000E+01- 

.700000E+01 
.100000E+01 

.300000E+01- 
. 000000E-99 

.100000E+01 
.100000E+01 
.300000E+01- 

.ooooooE-99- 

.300000E-S-01- 

. 000000E-99 
.000000E-99- 
300000E+01 
000000E-99 
.500000E+01- 



.300000E+01-, 

.280000E+02-, 

.ooooooE-99 
J90000E+02 

.400000E+01- 
. 120000E+02 
J00000E+01 

.360000E+02 

.400000E+01 

.420000E+02- 

.200000E+01 
.500000E+01 
.900000E+01 

.180000E+02 
.OOOOOOE-99 
. 190000E+02- 



,400000E+01 , 
J40000E+02-, 

, 000000E-99-. 
,95ooooe+oi . 
.500000E+01 , 

.600000E+01 , 
.100000E+01 . 
J80000E+02 , 

.000000E-99 . 

.210000E+02-. 
.200000E+01 . 
•250000E+01 , 

, 000000E-99 . 

.900000E+01 . 

, 000000E-99 , 
.950000E+01-. 



500000E+01- 

840000E+02 

200000E+01 

570000E+02 
000000E-99 
360000E+02 

100000E+01 

1 08000E+03 

100000E+01 
126000E+03 

000000E-99 
1 50000E+02 

700000E+01 

540000E+02 
000 00 oe -99 
570000E+02 



.600000E+01 

.600000E+02 
.000000E-99 
.300000E+02 
.500000E+01 
.900000E+02 

.100000E+01 

.450000E+02 

.200000E+01 

.135000E+03 

. 100000E+01 

.675000E+02 
.000000E-99 
.990000E+02 
.OOOOOOE-99 
.000000E-99 
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Problem No. 4 

13 18 

.OOOOOOE-99 .350000E+02 .100000E+02 .400000E+02 .100000E+02 .3OOO00E+02 

.150000E+03 .115000E+03 .4oooooe+o2 .410000E+02 .530000E+02 .550000E+02 

.67OOOOE+O2 .lUlOOOE+03 .170000E+02 .200000E+02 .380000E+02 .OOOOOOE-99 

-.350000E+02 .000000E-99 .nooooE+02 .150000E+02 .170000E+02 .130000E+02 

.120000E+02 .610000E+02 .640000E+02 .630000E+02 .720000E+02 .750000E+02 
.470000E+02 .480000E+02 .490000E+02 .520000E+02 .150000E+02-.700000E+02 
-.100000E+02-.1 10000E+02 .OOOOOOE-99 .100000E+01 .200000E+01 .300000E+01 

.4oooooe+oi .500000E+01 ..6oooooe+oi .700000E+01 .800000E+01 .900000E+01 

.100000E+02 .110000E+02 .120000E+02 .I3OOOOE+O2 . 140000E+02-.200000E+02 
-.400000E+02-.150000E+02-.100000E+01 .OOOOOOE-99 .110000E+02 . 1 50000E+02 
.320000E+02 .2l0000E+02 .170000E+02 .110000E+02 .150000E+02 . 160000E+02 
.170000E+02 .180000E+02 .190000E+02 .200000E+02 .2 1 OOOOE+02-.800000E+02 
-. 100000E+02-.170000E+02-.200000E+01-.1 10000E+02 .OOOOOOE-99 .I5OOOOE+O3 

.149000E+03 .148000E+03 .147000E+03 .146000E+03 .145000E+03 .144000E+03 
.143000E+03 .142000E+03 .141000E+03 .140000E+03 .139000E+03-.200000E+02 
-. 300000E+02-. i 30000E+02-. 300000E+0 1 -. 1 50000E+02-. 1 50000E+03 .000000E-99 

-. 200000E+0 1 -. 300000E+0 1 -. 400000E+0 1 -. 500000E+0 1 600000E+0 1 -. 700000E+0 1 
800000E+0 1 -. 900000E+0 1 -. 1 00000E+02-. 1 1 0000E+02-. 1 20000E+02-. 600000E+02 

-. 1 5OOOOE+03-. 12000OE+O2-.4OOO0OE+01-.3200O0E+02-.1490O0E+O3 .200000E+01 
.OOOOOOE-99 .210000E+02 .220000E+02 .230000E+02 .240000E+02 .250OO0E+O2 
.260000E+02 .270000E+02 .280000E+02 .290000E+02 .3OOOOOE+O2-.3OOOO0E+O3 

-. 1 I5OOOE+O3-.6IOOOOE+O2-.5OOOOOE+OI-.23OOOOE+O2-.1H80OOE+O3 .300000E+01 

-.210000E+02 .OOOOOOE-99 .310000E+02 .35OOOOE+O2 .39OOOOE+O2 .430000E+02 

.470000E+02 .147000E+03 .148000E+03 .149000E+03 .150000E+03-.230000E+03 

-. 400000E+02-. 640000E+02-. 600000E+0 1 -. 1 70000E+02-. 1 47OOOE+O3 . 400000E+0 1 
-.220000E+02-.310000E+02 .OOOOOOE-99 .160000E+03 .25OOOOE+03 .349OOOE+O3 
.47IOOOE+O3 .582OOOE+O3 .65IOOOE+O3 .78IOOOE+O3 ,8TiOOOE+03-.800000E+02 
-.410000E+02-.630000E+02-.700000E+01-.1 10OOOE+O2-.146OOOE+O3 .^OOOOOE+OI 

-.230000E+02-.350000E+02-. i6ooooe+03 . 000000E-99-. 300000E+02-.400000E+02 
-.500000E+02-.600000E+02-.700000E+02-.800000E+02-.900000E+02-.820000E+02 

-.53O00OE+O2-.720O00E+O2-.8O0000E+01-. 150OO0E+O2-.1450O0E+O3 .600000E+01 

-.240000E+02-.390000E+02— .250000E+03 .300000E+02 .OOOOOOE-99-. 130000E+02 

-. 1 50000E+02-. 140000E+02-. 160000E+02-. 1 50000E+02-. 170000E+02-. 106000E+03 
-.55OOO0E+O2-.75OOO0E+O2-.9OOOO0E+O1-. 160000E+02-. 144000E+03 .700000E+01 
-.250000E+02-.430000E+02-.349000E+03 .400000E+02 .130000E+02 .OOOOOOE-99 
.in000E+03 .121000E+03 . 133OOOE+O3 ..144000E+03 .I55OOOE+O3-.I 10000E+03 
-.67OOOOE+O2-.47OOOOE+O2-. 100000E+02-. 17OO00E+02-. I43OOOE+03 .800000E+0T 
-.260000E+02-.470000E+02-.471000E+03 .500000E+02 . 1 50000E+02-. 1 11000E+03 

.000000E-99 .nooooE+02 .900000E+01 .700000E+01 .500000E+01-.134000E+03 
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Problem No. 5 



13 18 

.000000E-99 
. 1 ^ooooe+O:! 
.670000E+02 
-.3<500ooe+o2 
. 120000E+02 
.470000E+02 

-. 100000E+Q2- 
.400000E+01 
. 100000E+02 
-.400000E+02- 
.320000E-K)2 
. 1 70000E+02 
-. 100000E+02- 
. 1 49000E+03 
.I43OOOE+03 
-.300000E+02- 
-.200000E+01- 
-.800000E+01- 
-. 150000E+03- 
.OOOOOOE-99 
.260000E+02 

-.1 15000E+03- 

-.210000E+02 
.470000E+02 
-.400000E+02- 
-.220000E+02- 

.47IOOOE+03 
-.410000E-J-02- 
-.230000E+02- 
-.500000E+02- 
-.53OOOOE-I-O2- 
-.240000E+02- 
-.1<50000E-i-02- 
-,5 I 50000E+02- 
-.250000E+02- 
.1 ' 1000E+03 
.670000E+02- 
-.260000E-I-02- 
.OOOOOOE-99 



.3 c )OOOOE+02 
.1 ^OOOE-K^ 
.14 1OOOE+O3 

.000000E-99 
.610000E+02 

.480000E+02 
.1 '0000E+02 
.500000E+01 
.1 10000E+02 
'.150000E+02- 
.230000E-f-02 
.180000E+02 
.170000E+02- 
. 1 48000E+03 
. 1 42000E+03 
..130000E+02- 
-.300000E+01- 
>.900000E+01- 
.120000E+02- 
.210000E+02 
.270000E+02 
.610000E+02- 

.000000E-99 
.147000E+03 

..640O0OE+02- 

■.310000E+02 
.582000E+03 
..630000E+02- 

..3 c 50000E-i-02- 
•.500000E+02- 
..720000E+02- 
■.390000E+02- 

• .140000E+02- 
..7^0000E+02- 
■.430000E+02- 
.121000E+03 
..470000E-I-02- 

• .470000E+02- 
.1 '0000E+02 



.100000E+02 
.400000E+02 
.I7OOOOE-I-O2 
. 1 10000E+02 
.640000E+02 
.490000E+02 
.OOOOOOE-99 
.60OOOOE+OI 
. 1 20000E+02 
.100000E+01 
. 17O000E+02 
.190000E+02 
.200000E+01- 
.I47OOOE+O3 
.141000E+03 
.3OOOOOE+OI- 
.400000E+01- 
.100000E+02- 
.400000E+01- 
.220000E+02 
.280000E+02 
.^OOOOOE+OI- 
.310O00E+O2 
.148000E+03 
.600000E+01- 

. 000000E-99 

^lOOOE-s-O^ 
.7OOOOOE+OI- 

. 160000E+03 

■.700000E+02- 
■■.OOOOOOE+01- 

..250000E+03 

.. 160000E+02- 
.900000E+01- 
..349000E+03 
.133000E+03 
..100000E+02- 
..471O0OE+03 
.9O0OO0E+01 



.400000E+02 
.410000E+02 
.200000E+02 
.1 t 50000E+02 
.630000E+02 
.520000E+02 
.100000E+01 
.7OOOOOE+OI 
. 1 30000E-I-02 
.OOOOOOE-99 
.1 ':0000E+02 
. 200000E+02 
. 1 0000E+02 
. 146000E-f-03 
. 1 l+OOOOE+03 
. 1 '50000E+02- 
.500000E+01- 
. 1 0000E+02- 
.320000E+02- 
.230000E+02 
.290000E+02 
.230000E+02- 
.3 t 5OO00E+O2 

.149000E+03 
. 170000E+02- 

.160000E+03 

.78ioooe-:-03 

. 1 1 0000E+02- 

.OOOOOOE-99- 

.800000E-!-02- 
. 1 ^OOOOE+02- 

. 300000E+02 

. 150000E-!-02- 
. 160000E+02- 
.400000E+02 
.144000E-I-03 
.. 170000E+02- 
.500000E+02 
.7OOOOOE+OI 



. 1 00000E-I-02 
. c 530000E+02 
.380OOOE+02 
. 170000E+02 
.720000E+02 
.1 ^OOOOE-i-02- 
.200000E+01 
.800000E-K)1 
.140000E+02- 
. 1 'OO0OE+O2 
. 1 c )OO0OE-!-O2 
.210000E+02- 
. OOOOOOE-99 
. A5OOOE+O3 
. 1 39OOOE+O3- 
. 1 ^OOOOE+O^ 
.600000E+01- 
■. 120000E+02- 
. 1H9OOOE+O3 
.240000E+02 
.300000E+02- 
. 148000E+03 
.390O0OE+02 
.1W00E+03- 

.147000E+03 
.2SO000E+03 
.81 IOOOE+03- 
. 146000E+03 
.300000E+02- 
.900000E+02- 
. 145000E+03 
.OOOOOOE-99- 
.. 170000E+02- 

.144000E+03 

. 1 30000E+02 
. 1 55000E+03- 
..i53oooe-k)3 
.1 50000E+02- 
. 500000E-K) 1 - 



.300000E-S-02 

.5 c )0000E+02 

.OOOOOOE-99 
. 1 30000E-5-02 
.750000E+02 

.700000E-S-02 

.300000E+01 
.900000E+01 

.200000E+02 

. 1 50000E+02 
.160000E+02 
.800000E+02 

.T : )OOOOE+03 

. 1 44000E+03 

.200000E+02 
.OOOOOOE-99 
■.700000E+01 
■.bOOOOOE+02 
.200000E+01 
,2^0000E+02 
.30OOOOE+03 
.300000E-S-01 
.430000E+02 
.230000E+03 
.400000E+01 
.3U9OOOE+03 
.800000E+02 
.S00000E-MD1 
.400000E+02 
^.820000E+02 
.600000E+01 
. 130000E+02 
•.106000E+O3 
.700000E+01 
.OOOOOOE-99 
..1 OOOOE+03 
.80O0O0E+01 
..1 1 1000E-5-03 
..134000E+03 
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Problem No. 6 

17 18 

.OOOOOOE-99 .3^0000E+02 .100000E+02 .^00000E+02 . 10000OE+O2 .3O00O0E+02 

.ISOOOOE+O^ .1 c 5000E+03 .bOOOOOE+02 .lHOOOOE+02 .53OOOOE+O2 .550OOOE+02 

.670OOOE+O2 .T^lOOOE+03 .170000E+02 .200000E+02 .380000E+02 .OOOOOOE-99 

-.350000E+02 .OOOOOOE-99 ,l:0000E+02 . 1 500O0E+02 .170O0OE+O2 .130000E+02 

.120000E+02 .610000E+02 ,6if0000E+02 .630OOOE+O2 .720000E+02 .750000E+02 

.U7OOOOE+O2 .480000E+02 .490000E+02 .520000E+02 .I5OOOOE+O2-.7OOOOOE+O2 

-.100000E+02-.1 10000E+02 .OOOOOOE-99 .100000E+01 .200000E+01 .3OOOOOE+OI 

.U00000E+01 .500000E+01 .600000E+01 .70ooooe+oi .800000E+01 .900000E+01 

.100000E+02 .110000E+02 .120000E+02 .130000E+02 . 140000E+02^.200000E+02 

-.400000E+02-.150000E+02-.100000E+01 .OOOOOOE-99 .110000E+02 .15000OE+02 

.320000E+02 .230000E+02 .170000E+02 . 110000E+02 .150000E+02 .1600O0E+O2 

. 1 70000E+02 .l80000E+02 .190000E+02 .200000E+02 .21 00OOE+O2-.80O0OOE+O2 

-.100000E+02-.170000E+02-.200000E+01-.1 10000E+02 .OOOOOOE-99 .I5OOOOE+O3 

. 1 49000E+03 . 1 48000E+03 . 1 47000E+03 . 1 46000E+03 . 1 45OOOE+O3 . 1 44000E+03 

.143000E+03 .142000E+03 .141000E+03 .140000E+03 .139000E+03-.200000E+02 

300000E+02-. 1 30000E+02-.300000E+01 -. 1 50000E+02-. 1 50000E+03 . OOOOOOE-99 

-.200000E+0 1 -.300000E+01 -.400000E+01 -. 500000E+01 -.600000E+01 -.700000E+01 
-. 800000E+0 1 -. 900000E+0 1 -. 1 00000E+02-. 1 1 00OOE+02-. 1 20000E+02-.600000E+02 
-. 1 5OOOOE+O3-. 120000E+02-.400000E+01-.320000E+02-. 149000E+O3 .200000E+01 
.OOOOOOE-99 .210000E+02 .220000E+02 .230000E+02 .240000E+02 .250000E+02 
.260000E+02 .270000E+02 .280000E+02 .2900O0E+O2 .3OO00OE+O2-.3OOOOOE+O3 
-.1 15OO0E+03-.61O00OE+O2-.500000E+01-.230000E+02-.148000E+03 .300000E+01 
-.210000E+02 .OOOOOOE-99 .310000E+02 .350000E+02 .390000E+02 .430000E+02 
.470000E+02 .147000E+03 .148000E+03 .149000E+03 .I50OOOE+O3-.23OOOOE+O3 
-.400000E+02-.640000E+02-.600000E+01-. 170000E+02-. I47OOOE+O3 .400000E+01 

-.220000E+02-.310000E+02 .000000E-99 .160000E+03 .250000E+03 .349000E+03 

.47IOOOE+03 .5820OOE+O3 .651000E+03 .78IOOOE+03 .81 IOOOE+O3-.80OOOOE+O2 
-.410000E+02-.630000E+02-.700000E+01-.1 10000E+O2-. 1460O0E+O3 .5OOOOOE+OI 

-.230000E4-02-.350000E+02-.I60000E+03 . 000000E-99-. 300000E+02-.400000E+02 
-.500000E+02-.600000E+02-.700000E+02-.800000E+02-.9000COE+02-.820000E+02 

-.530000E+02-.720000E+02-.500000E+01-.150000E+02-. 145000E+03' .60000OE+01 
-.240000E+02-.390000E+02-.250000E+03 .300000E+02 .OOOOOOF 99-.I3OOOOE+O2 
-. 1 50000E+02-. 1 40000E+02-. 1 60000E+02-. 1 5000OE+O2-. 1 700OOE+O2-. 1 060OOE+O3 
-.550000E+02-.750000E+02-.900000E+01-. 16O000E+O2-.1440OOE+03 .7OOOOOE+OI 
-.250000E+02-.430000E+02-.349000E+03 .400000E+02 .130000E+02 .OOOOOOE-99 

.1-1000E+03 .121000E+03 .133000E+03 .144000E+03 . 155000E+03-.1 100OOE+03 
-.670000E+02-.470000E+02-. 100000E+02-. 170000E+02-. I43OOOE+O3- ^800000E+O1 
-.26O000E+02-.4700O0E+02-.47100OE+O3 .5000O0E+O2 . 1 50000E+02-. 1 IIOOOE+03 

.OOOOOOE-99 .T0000E+02 .9OOOOOE+OI .700000E+01 . 500000E+01 -. 1 34000E+03 
141000E+03-.480000E+02-. 1 10000E+02-. l80000E+02-. 142000E+03 .900000E+01 
-. 270000E+02-. 1 47000E+03-. 5820OOE+O3 .600000E+02 . 1 40000E+02-. 1 21 OOOE+03 

-.nooooE+02 .000000E-99 .150000E+02 .200000E+02 .300000E+02-.282000E+03 

-. 1 70000E+02-.490000E+02-. 1 20000E+02-. 1 90000E+02-. 141 000E+03 . 1 000OOE+02 
-.280000E+02-. 1 480OOE+O3-.65I OOOE+03 .700000E+02 . 1 60000E+02-. 1 33OOOE+O3 
-.900000E+01-.150000E+02 .OOOOOOE-99 .IOOOOOE+01 .20OOOOE+O2-.34OOOOE+O2 
-.20OOOOE+O2-.52O0OOE+O2-. 130000E+02-.200000E+02-. 140000E+03 . 1 1000OE+O2 
-.290000E+02-. 149000E+03-.781000E+03 .800000E+02 . 1 50000E+02-. 144000E+03 
-.700000E+01-.200000E+02-. 100000E+01 .OOOOOOE-99 . 100000E+02-.400000E+02 
-,38o0O0E-i-O2-.15O0OOE+O2-.14O0OOE+O2-.21O000E+02-.13900OE+O3 .120000E+02 
-.3OO0OOE+O2-. I5OOOOE+O3-.8I OOOE-J-03 .900000E+02 . 170000E+02-. I55OOOE+O3 
-. 500000E+0 1 -. 300000E+02-. 200000E+02-. 1 00000E+02 . OOOOOOE-99-. 7bOOOOE+02 




Problem No. 7 

17 18 

.OOOOOOE-99 
. 1 5OOOOE+O3 
.67OOOOE+02 
■.350000E+02 
. 120000E+02 
.470000E+02 

• . 100000E+02- 
•400000E+01 
. 100000E+02 

• .400000E+02- 
.320000E+02 
.170000E+02 

-. 100000E+02- 
. 149OO0E+O3 
. 143000E+03 
-.3OOOO0E-fO2- 
•.200000E+01- 
■.800000E+01- 

■ .150000E+03- 

. 000000E-99 
.260000E+02 
..115000E+03- 

• .210000E+02 
.470000E+02 
..400000E+02- 

• .220000E+02- 
.U7IOOOE+03 

..410000E+02- 
-.230000E+02- 
■.500000E+02- 
-.530000E+02- 
-.240000E+02- 
■.150000E+02- 
-.550000E+02- 
-.250000E+02- 

.111 000E+03 
-.670000E+02- 
-.260000E+02- 

.OOOOOOE-99 
-.141000E+03- 
-.270000E+02- 
-.1 10000E+02 
-. 170000E+02- 
-.280000E+02- 
-.900000E+01- 
-.200000E+02- 
-.290000E+02- 
-.700000E+01- 

.380000E-I-02- 
..300000E+02- 

■ .500000E+01- 



• 3 I 50000E+02 
. 1 1 5000E+03 
.lHl000E+03 
.OOOOOOE-99 
.610000E+02 
.480000E+02 
. 1 10000E+02 
.500000E+01 
. 1 10000E+02 
150000E+02- 
.230000E+02 
. 180000E+02 
. 170000E+02- 
. 148000E-:-03 
. 142000E+03 
130000E+02- 
.300000E+01- 
.900000E+01- 
. 120000E+02- 
.210000E+02 
.270000E+02 
.610000E+02- 
.OOOOOOE-99 
.147000E4-03 
.640000E+02- 
.310000E+02 
.5820OOE+O3 
.630000E+02- 
.350000E+02- 
.600000E+02- 
.720000E-J-02- 
.390000E+02- 
.140000E+02- 
.750000E+02- 
.430000E+02- 
.121000E+03 
.470000E+02- 
.470000E+02- 
. 1 1 0000E+02 
.480000E+02- 
.I47OOOE+O3- 

. 000000E-99 

.490000E+02- 
. I480OOE+O3- 
. 1 50000E+02 
.520000E+02- 

.149000E+03. 
.200000E+02- 
. 15000 0E+02- 

. 1 c )0000E+03- 

.300000E+02- 



. 100000E+02 
.400000E+02 
. 170000E+02 
. 1 10000E+02 
.640000E-f02 
.^90OO0E-fO2 
.OOOOOOE-99 
.600000E+C1 
.120000E+02 
. 100000E+0 1 
. 170000E+02 
. 190000E+02 
.200000E+01- 
. 147000E+03 

.141000E+03 
.300000E+01- 

.400000EH-01- 

.100000E+C2- 

.400000E+01- 

.220000E+02 

.280000E+02 

.500000E+01- 

.310000E+02 

.148000E+03 

.600000E-:-01- 

.OOOOOOE-99 

.65IOOOE+03 

.700000E+01- 

.160000E-S-03 
.700000E+02- 
.800000E+01- 
.250000E+03 

160000E+02- 
.900000E+01- 
■.349000E+03 
.133O00E+03 
.100000E+02- 
.47IOOOE+O3 
.900000E+01 
. 1 10000E+02- 
..582000E+03 
. 150000E+02 
.. 120000E+02- 
..65IOOOE+O3 

.■ooooooE-99 
..130000E+02- 
..781000E+03 
..100000E+01 

..140000E+02- 
..81 1000E+03 
.200000E-:-02- 



.400000E+02 
,Ul0000E-f02 
.200000E+02 
. 1 50000E+02 
.D30000E+02 
•520000E+02 
, 100000E+01 
.700000E+01 
. 130000E+02 

.ooooooe-99 
, 1 -0000E+02 

.200000E+02 
. 1 " 0000E+02 
, 14oOOOE+03 
. 14OO00E+O3 
. 150000E+02- 
,500000E-!-01- 
, 1 10000E+02- 
.320000E+02- 
.230OO0E+O2 
.290000E+02 
.230000E+02- 
.350000E+02 
, 159000E+03 
, 170000E+02- 
.160000E+03 

•78ioooe+03 

. 1 10000E+02- 

.ooooooe-99- 

.800000E+02- 

. 1 50000E+02- 
.300000E+02 
, 150000E+02- 
, 160000E+02- 

.400000E+02 

, 1 44000E+03 
. 170000E+02- 
, 500000E+02 
.700000E+01 

,l80000E+02- 

.600000E+02 

.200000E+02 

. 190000E+02- 

.700000E+02 

, 100000E+01 

.200000E+02- 

.800000E+02 

.OOOOOOE-99 

.210000E+02- 

, 900000E+02 

.100000E+02 



100000E+02 

530000E+02 

38OO00E+O2 

170000E+02 

720000E+02 

150000E+02- 

200000E-S-01 

800000E+01 

140000E+02- 

1 • OOOOE+02 

1 50000E-I-02 

210000E+02- 

OOOOOOE-99 

I45OOOE+O3 

139000E+03- 

1 5OOOOE+O3 

600000E+01- 

120000E+02- 

149000E+03 

240000E+02 
300000E+02- 

148000E+03 

390000E+02 

150000E+03- 

147000E+03 

250000E+03 

8i 1000E+03- 

146000E+03 

300000E+02- 

90000 0E+02 - 

145000E+03 

000000E-99- 

170000E+02- 

i^UoooE+03 

1 30000E+02 

1 55000E+03- 

143000E+03 

150000E+02- 

500000E+01- 

142000E+03 

1U0000E+02- 

300000E+02- 

141000E+03 

160000E+02- 

200000E+02- 

140000E+03 
150000E+02- 

100000E+02- 
139000E+03 
170000E+02- 
OOOOOOE-99- 
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.300000E+02 
.550000E+02 
.OOOOOOE-99 
. 130000E+02 
.750000E+02 
700000E+02 
300000E+01 
. 90000 0E+01 
.200000E+02 
.150000E+02 
. 1 60000E+02 
.800000E+02 
. 1 5OOOOE+O3 
. 144000E+03 
.200000E+02 

.000000E-99 

.700000E+01 
.600O00E+02 
.200000E+01 
.250000E+02 
.3OOOOOE+03 
.300000E+01 
.430000E+02 
.230000E+03 
.400000E+01 
.349OOOE+O3 
00000E+02 
.S00000E+01 
.400000E+02 
.820000E+02 
.600000E+01 
.130000E+02 
. 106000E+03 
• 700000E-f01 
.OOOOOOE-99 
.1 10000E+03 
.800000E+01 
. 1 1 1000E+03 
.134000E+03 
.900000E+01 
. 121000E+03 
.282O00E+O3 
. 1 00000E+02 
. 133000E+03 
.340000E+02 
. 1 ■ 0000E+02 
.144000E+03 
.400000E+02 
. 120000E+02 
. 1 55OOOE+O3 
.760000E+02 



a 



Page 41 

Comments on the Typewriter Log for Sample Problems 



Pages 42-50 include the typewriter log obtained when the program 
was used to solve the sample problems. The following should be 
noted: 

(a) A tolerable zero value was entered for each one of the 
problems . 

(b) Problems No. 1, 2, 4, and 6 each has an infinite number 
of solutions; SW3 was set to the ON position in order to 
obtain the general solution for each one of them. 

(c) The tolerable Z value for Problem No. 2 was typed incor- 
rectly; SW4 was set to the ON position, and the Z value 
was re-entered in the proper way. 

(d) The tolerable Z value for Problem No. 4 was typed incor- 
rectly and SW4 was left in the OFF position. The message 
ERROR F7 was typed and the program continued. According 
to the OPERATING INSTRUCTIONS, the program was stopped and 
re-initialized. 

(e) Problems No. 5, 6, and 7 each required readjustment of the 
tolerable Z value. This was done according to instructions 



Typewriter Log for Sample Problems 



160001000000 
LOAD DATA 

ENTER TOLERABLE Z IN E\k.8 FORMAT SPEC. 
+.10000000E-05 

SOLUTION FOR VECTOR 5 

INF SOLS SW3 ON FOR GEN S0L-SW3 OFF TO CONT. 
X 3 

-.10000000E+01 
.10000000E+01 

X k 
.00000000E-99 
.10000000E+01 



.30000000E+01 
-.20000000E+01 

ENTER TOLERABLE Z IN E 1 if . 8 FORMAT SPEC. 

+.10000000E-OJ 

+ .10000000E-0if 

SOLUTION FOR VECTOR 6 

INF SOLS SW3 ON FOR GEN S0L-SW3 OFF TO CONT. 
X 5 

-.87M+6120E+01 
-. 17^72361 E+01 
.2211 1665E+01 
-.15055279E+01 
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-.OO000O00E-99 
-.O0OO0OO0E-99 
-.00000000E-99 
. 00000000E-99 

ENTER TOLERABLE Z IN E1*f.8 FORMAT SPEC. 
+.10000000E-03 

SOLUTION FOR VECTOR 9 
. 9999997 OE-00 
.19999998E+01 
.29999999E+OI 
.if0000002E+01 
. i+999999 0E+01 
.6000001 OE+01 
. 69999996 E+01 
.79999996E+01 

SOLUTION FOR VECTOR 10 
. 4999997 OE-00 
. 1 0000003E+01 
.15000001E+01 
.20000000E+01 
.24999999E+01 
.3 0000005 E+01 
.34999996E+01- 
•39999995E+01 

SOLUTION FOR VECTOR 11 
.299999&7E+01 
.5999999^+01 
.89999980E+01 
.11999999E+02 
.1 4999995 E+02 
.1800O0OCE+02 
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•2 0999998E+02 
. 23999996 E+02 

SOLUTION FOR VECTOR 12 

•17373927E+02 

.46M7256E+01 
-. 1 1788599E+02 

.23638003E+01 

.11300876E+02 

. 9985957 OE-OO 

.125685^5E+02 

.75^1 1273E+01 

ENTER TOLERABLE Z IN E14.8 FORMAT SPEC. 

+.10000000E-00 

ERROR F7 

4908000 

ENTER TOLERABLE Z IN E14.8 FORMAT SPEC. 
+.10000000E-03 

SOLUTION FOR VECTOR ]k 

SYSTEM INCONSISTENT 

SOLUTION FOR VECTOR 15 

SYSTEM INCONSISTENT 

SOLUTION FOR VECTOR 16 

SYSTEM INCONSISTENT 

SOLUTION FOR VECTOR 17 



SYSTEM INCONSISTENT 



SOLUTION FOR VECTOR 18 



INF SOLS SW3 ON FOR GEN SOL-SW3 OFF TO CONT. 
X 13 

-.89912030E-00 
-.11049142E+01 

•57970323E+02 
-.9^171 i t20E+01 

.1 kl 4684 OE-OO 

. 17722926E+01 

. 752997 10E+01 
-.1 23M+7 10E+02 

. 1 1761880E+01 

. 6271787 OE+01 
-. 356827 07E+01 
-.20l663ifOE+01 



.20000000E+01 
-.OOOOOOOOE-99 
-.OOOOOOOOE-99 
-.OOOOOOOOE-99 
-, OOOOOOOOE-99 
-.OOOOOOOOE-99 
-.OOOOOOOOE-99 
-.OOOOOOOOE-99 
-.OOOOOOOOE-99 
-.OOOOOOOOE-99 
-.OOOOOOOOE-99 

.OOOOOOOOE-99 

ENTER TOLERABLE Z IN E14.8 FORMAT SPEC. 
+.10O00O00E-03 

RELOAD DATA I NCR. Z BY FACTOR OF 10 
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ENTER TOLERABLE Z IN E14.8 FORMAT SPEC, 
+.10000000E-02 

SOLUTION FOR VECTOR 14 

-.483861 03E+01 

. J+1 42566 1E+03 
-.55^8 039 1E+02 

. 1427631 1E-00 

.15^99623E+02 

.56783579E+02 
-.996295 10E+02 

.1031 0057 E+02 

.47450099E+02 
-.12369882E+02 
-.26474642E+02 

.91686745E+01 

SOLUTION FOR VECTOR 15 

-.77705800E+01 

-.49968307E+01 

. 37^0047 OE+03 
-.51928873E+02 

.33781753E-00 

.15138997E+02 

.49755229E+02 
-.88164190E+02 

.884691 10E+01 

.39042339E+02 
-.89736869E+01 
-.23229999E+02 

.81447152E+01 

SOLUTION FOR VECTOR 16 
-.7037 1 *3^0E+01 
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-.44552765E+01 

.33750718E+03 
-.46798547E+02 

.29713861E-00 

.13666631 E+02 

.45032 153 E+02 
-.79983800E+02 

.79584935E+01 

.34742553E+O2 
-.74102383E+01 
-.20863395E+02 

.74716143E+01 

SOLUTION FOR VECTOR 17 

-.687297^0E+01 

-.50068592E+01 

.3^87341 OE+03 
-.50506716E+O2 

. 45286507 E-00 

.121750O5E+02 

.47732087E+02 
-.82814740E+02 

.8401 823 OE+01 

.38005589E+02 
-.9595 0557 E+01 
-.2 1404485 E+02 

.79367011 E+01 

SOLUTION FOR VECTOR 18 

.96000000E-05 
-.24577692E+01 

.12894879E+03 
-.20947452E+02 

.32802864E-00 

.39422583E+01 

.1 67496 22E+02 
-.27459^73E+02 



.26163002E+01 

. 139509 05E+02 
-.79372437EH-01 
-.448577 05 E+01 

.22243907E+01 

ENTER TOLERABLE Z IN E14.8 FORMAT SPEC. 
+. 10000000E-03 

RELOAD DATA I NCR. Z BY FACTOR OF 10 

ENTER TOLERABLE Z IN E14.8 FORMAT SPEC. 
+.10000000E-02 

SOLUTION FOR VECTOR 18 

INF SOLS SW3 ON FOR GEN S0L-SW3 OFF TO CONT. 
X 17 

-.42874995E-00 

.23976300E+01 
-.28644234E+02 

.97720987E+01 
-.61967767E-00 
-.17527315E+01 
-.82128261E-00 

. 2129996 0E-01 
-.25086718E-00 
-. 18757 052E+01 

.35436282E+01 

.19 5 03 036E+01 

.4330901 1E-00 
-.18223183E+01 

.10547266E+02 
-.10317607E+02 



.20000000E+01 



-.00000000E-99 
-.00000000E-99 
-.00000000E-99 
-.OOO00O0OE-99 
-.00000000E-99 
-.O0O0OO0OE-99 
-.OOOOOOOOE-99 
-.00000000E-99 
-.00000000E-99 
-.OOOOOOOOE-99 
-.OOOOOOOOE-99 
-.OOOOOOOOE-99 
-. OOOOOOOOE-99 
-.OOOOOOOOE-99 
.OOOOOOOOE-99 

ENTER TOLERABLE Z IN E14.8 FORMAT SPEC. 
+.10000000E-03 

RELOAD DATA I NCR. Z BY FACTOR OF 10 

ENTER TOLERABLE Z IN E14.8 FORMAT SPEC. 
+.10000000E-02 

RELOAD DATA I NCR. Z BY FACTOR OF 10 

ENTER TOLERABLE Z IN E 14.8 FORMAT SPEC. 
+.10000000E-01 

SOLUTION FOR VECTOR 18 

.91800000E-04 

.11 18371 9E+02 
-.13361065E+03 

.45581915E+02 
-.28904901 E+01 
-.81754060E+01 
-.38309755E+01 

.99337500E-01 
-.1 1702277E+01 
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-.87^971 3 0E+01 
.16529699E+02 
.90972961E+01 
.20201331E+01 

-.85001620E+01 
.if9197987E+02 

-.481 269^5 E+02 
.466463 03E+01 




Detailed Flow Chart 
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d i me ns i on 
a(T7,i8),pvr(78),x(T7), 

tr(78) 




1=1 








J=1 



Mote; 

All exits and en- 
trances at the bottom 
of a page are concerned 
only with entrances sad 
exits, respectively, at 
the top of the following 
page. 



(k{ I,J),a( I,J+1),a( I,J+2),a( I, J+3),a( l,J+4), 
a( l,J+5) 

;ormat(E12.6,E12.6,E12.6,E12.6,E12.6,E12.6) 
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Program Listing 

38000 C SOLUTION OF A SYSTEM OF N LINEAR EQUATIONS WITH N UNKNOWNS 
38000 C BY THE METHOD OF ELEMENTARY OPERATIONS. 
38000 C ENGINEER *HEBERTO PACHON 



08000 




DIMENSION A(17,18),PVR(18),X(17),TR(18) 


38000 


1 


READ 50,N,M 


38036 




DO 2 1=1 ,N 


38048 




DO 2 J=1,M,6 


38060 


2 


READ 51 ,A( 1 , J) ,A( 1 , J+1 ) ,A( 1 ,J+2) ,A( 1 , J+3) ,A( 1 , J+4) ,A( 1 ,J+5) 


38576 




PRINT 52 


38600 




ACCEPT 53, Z 


38624 




DO 3 1=1 ,N 


38636 


3 


X(l)=0. 


38720 




CV=0. 


38744 




KC=0 


38768 




MC=0 


38792 




NR=0 


38816 


4 


MC=MC+1 


38852 




NR=NR+1 


38888 


5 


IF(MC-N)6,6,22 


38956 


6 


K=NR 


38980 


7 


IF(K-N)8, 10,10' 


39048 


8 


KK=K+1 


39084 




DO 9 l=KK,N 


39096 




IF(ABS(A(K,MC))-ABS(A(I ,MC))) 12,9,9 


39320 


9 


CONTINUE 


39356 


10 


IF(ABS(A(K,MC))-Z)11,11,13 


39484 


11 


KC=KC+1 


39520 




X(MC)=1. 


39568 




MC=MC+1 



39604 


GO TO 5 




39612 12 


K=l 




39636 


GO TO 7 




39644 13 


DO 15 J=MC,M 




39656 


PVR(J)=A(K,J) 




39764 


A(K,J)=A(NR,J) 




39908 


A(NR,J)=PVR(J) 




T0016 


IF(MC-1)15,14,15 




T0084 14 


TR(J)=PVR(J) 




T0156 15 


CONTINUE 




T0192 


PVR(NR)=A(NR,MC) 




T0300 


DO 16 J=MC,M 




T0312 16 


A(NR,J)=A(NR,J)/PVR(NR) 


T0528 


IF(MC-M) 18, 17,61 




T05 96 1 7 


IF(NR-N) 18,20,61 




T0664 18 


L=NR+1 




T0700 


DO 19 l=L,N 




T0712 


PVR(I)=A(I,MC) 




T0820 


DO 19 J=MC,M 




T0832 19 


A(l ,J)=A(I ,J)-PVR( 


)*A(NR,J) 


T1168 


IF(MC-1)20,4,20 




T1236 20 


L=NR-1 




T1272 


DO 21 1=1, L 




T1284 


PVR(I)=A(I ,HC) 




T1392 


DO 21 J=MC,M 




T1404 21 


A(l ,J)=A(I ,J)-PVR( 


)*A(NR,J) 


T1740 


GO TO 4 




T1748 22 


IF(KC)23,37,23 




T1804 23 


NR=N-KC+1 




T1852 


KR=N-KC 




T1888 


L=N+1 





T1924 


00 2k 1 = 1, KR 


Tl 936 2k 


PVR( 1 )=A( 1 ,L) 


T2080 


00 27 J=l ,N 


T2092 


IF(X(J))25,27,25 


T2172 25 


CV=CV+TR(J) 


T2232 


KK=J-1 


T2268 


NL=J 


T2292 


DO 26 K=1 ,KK 


T2304 


TR( NL)=TR( NL-1 ) 


T2376 26 


NL-NL-1 


T2M*8 


DO 27 1=1, KR 


T2460 


PVR(I)=PVR(!)-A(I,J) 


T2604 27 


CONTINUE 


T2676 


DO 28 1=1, KR 


T2688 


KK= 1 +KC 


T272^ 28 


CV=CV+TR(KK)*PVR( 1 ) 


T2856 


IF(ABS(TR(L)-CV)-Z)29,29,59 


T2972 29 


PRINT 5k,L 


T2996 


DO 30 l=NR,N 


T3008 


IF(ABS(A(I ,L))-Z)30,30,36 


T3136 30 


CONTINUE 


T3172 


PRINT 55 


T3196 


PAUSE 


T3208 


IF(SENSE SWITCH 3)31,35 


T3228 31 


DO 33 J=1,N 


T3240 


IF(X(J))32,33,32 


T3320 32 


PRINT 56, J 


T33M* 


DO 33 1=1, KR 


T3356 


A(I,J)=-A(I,J) 


T3512 


PRINT 53,A(I,J) 


T3596 


A(I,J)=-A(I,J) 
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T3752 


33 


CONTINUE 


T382** 




PRINT 57 


T3848 




DO 3k i=l,KR 


73860 


3k 


PRINT 53,A{i,L) 


T3980 


35 


L=L+1 


T^016 




IF(L-M)29,29,61 


lk08k 


36 


PRINT 58 


Ik] 08 




GO TO 35 


Ti+1 16 


37 


J=N+1 


T4152 




DO 38 l=1,N 


lk\6k 


38 


CV=CV+TR( 1 )*A( 1 ,J) 


T4332 




IF(ABS(TR(J)-CV)-Z)39,39,59 


lkkk8 


39 


PRINT 5k, J 


lkk72 




DO kO 1=1, N 


lkk8k 


ko 


PRINT 53, A(l, J) 


IkSok 




J=J+1 


Ik6k0 




IF(J-M)39,39,61 


T4708 


50 ■ 


FORMAT ( 13,13) 


T4736 


51 


F0RMAT(E12.6,E12.6,E12.6,E12.6,E12.6,E12.6) 


lkj8k 


52 


FORMAT(/39HENTER TOLERABLE Z IN E 1 ^ . 8 FORMAT SPEC.) 


lk892 


53 


F0RMAT(E14.8) 


IkS Ik 


Sk 


F0RMAT(/19HS0LUTI0N FOR VECTOR, I 3) 


T4986 


55 


FORMAT( A4HINF SOLS SW3 ON FOR GEN S0L-SW3 OFF TO CONT.) 


T5104 


56 


F0RMAT(/1HX, 1 3) 


T5140 


57 


F0RMAT(/) 


T5162 


58 


FORMAT(/ 19HSYSTEM INCONSISTENT) 


T5230 


59 


PRINT 60 


T525^ 


60 


FORMAT( 35HRELOAD DATA INCR. Z BY FACTOR OF 10) 


T5348 




PAUSE 


T5360 


61 


GO TO 1 


T5368 




END 
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T9999 


SIN 




T6139 


0052 


T9989 


SINF 




T6129 


5053 


T9979 


COS 




T6119 


0053 


T9969 


COSF 




T6109 


Z 


T9959 


ATAN 




T6099 


(3003 


T9949 


ATANF 




T6089 


5ooooooo55 


T9939 


EXP 




T6079 


cv 


T9929 


EXPF 




T6069 


KC 


T9919 


LOG 




T6059 


5oo5 


T9909 


LOGF 




T6049 


MC 


T9899 


SQRT 




T6039 


NR 


T9889 


SQRTF 




T6029 


5oo4 


T9879 


ABS 




T6019 


5ooT 


T9869 


ABSF 




T6009 


500 


T9859 


RND 




T5999 


(5005 


T9849 


RNDF 




T5989 


5oo6 


T9839 


A 


T6789 


T5979 


5022 


T6779 


PVR 


T6609 


T5969 


K 


T6599 


X 


T6439 


T5959 


T5007 


T6429 


TR 


T6259 


T5949 


5008 


T6249 


"5001 




T5939 


5010 


T6239 


■D050 




T5929 


KK 


T6229 


T5050 




T5919 


5009 


T6219 


N 




T5909 


501 


T6209 


M 




T5899 


5012 


T6199 


0002 




T5889 


5011 


T6189 


1 




T5879 


5013 


T6179 


J 




T5869 


T000000051 


T6169 


0051 




T5859 


5015 


T6159 


0051 




T5849 


5oi if 


T6149 


1)052 




T5839 


5016 



T5829 5018 


T5519 5057 


T5819 5017 


T5509 5034 


T5809 5061 


T5499 5058 


T5799 5020 


T5489 5058 


T5789 L 


T5479 5038 


T5779 5019 


T5469 5039 


T5769 5021 


T5459 5040 


T5759 5023 


T5449 5060 


T5749 5037 


T5439 5060 


T5739 KR 




T5729 5024 




T5719 5027 




T57 09 5025 




T5699 NL 




T5689 5026 




T5679 5028 




T5669 5029 




T5659 5059 




T5649 5054 




T5639 5054 




T5629 5030 




T5619 5036 




T5609 5055 




T5599 5055 




T5589 5031 




T5579 5035 




T5569 5033 




T5559 5032 




T5549 5056 




T5539 5056 




T5529 5057 
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Maintenance 



I . Core Layout 

The program listing and the symbol table show the storage addresses 
for the statements, variables, subroutines, and constants used in 
conjunction with the program. It is noted that 

(a) Program starts at location 08000 

(b) The five digit address on the left of each source statement 
is the starting address in core storage for instructions 
compiled for that statement. 

(c) The address 15368 associated with the END statement is 
the first available location after the last program 
instruction, excluding subroutines. 

(d) Sixteen ten digit fields are located in positions 19840 
through 19999. The alphabetic representation of the 
names of the eight relocatable subroutines, in the two 
forms allowed — one with and one without the terminal F — 
are stored in the 16 fields . 

(e) Unsubscripted variables used in the source program 
appear as 

Field address Variable name 

(f) Subscripted variables appear as 

Field address Variable name Field address 
first element last element 

(g) Floating point constants appear as 

Field address constant (eight digit mantissa and 
two digit exponent) 

(h) Fixed point constants appear as 

Field address constant (four digits right 
justified) 
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(i) Statement numbers appear as 

Field address Statement 
of compiled branch number 

Fields are addressed from their units position (high core position). 

The first 8000 positions include the multiply-add table?, and the 
arithmetic and input/output subroutines together with the work areas 
they require. 



II. Method of Incorporating Changes 

If the programmer wishes to use this program for a system of equations 
that do not require the maximum storage, he may economize core storage 
by reducing the dimensioned areas A(17,18), PVR(18), X(17), TR(18) 
appearing in the first source program statement. 

The following rules should be observed: 

(a) Determine the number N 

N = number of equations =» number of unknowns 

(b) Determine the number M 

M = number of equations plus number of constant 
vectors for which solutions are desired 



(c) Adjust the DIMENSION statement according to the table 
below: 



Value of M 


DIMENSION STATEMENT 


M £ 6 


DIMENSION A(N, 6), PVR(6), X(N), TR(6) 


6 < M s 12 


DIMENSION A(N, 12), PVR(12), X(N) , TI<L2) 


12 < M s 18 


DIMENSION A(N, 18). PVR(18), X(N) , TR(18) 
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This method of reducing core storage was used in conjunction with a 
problem which required the solution of a system of five equations 
with five unknowns and one constant vector. The specifications 
further required a repeated solution of the system for each of the 
coefficient matrices resulting from the change of a variable parameter 
affecting the matrix elements . 

To accomplish this problem, the DIMENSION statement was adjusted ac- 
cording to the above specifications , and the loop for reading the 
coefficient elements (DO 2 1= 1,N; DO 2 J = 1,M; 2 READ 51, A(I,J), 
etc., etc.,....) was replaced by a recursive routine which would 
generate the coefficient matrix elements, each time the variable 
parameter changed. The program would then evaluate the unknowns 
for each new system so generated. 



III. Hash Totals 

Two hash totals obtained by means of the CARD HASH TOTAL program 
written by Mr. William G. Weideman of the Marquette University 
Computer Center are listed below: 

(1) Hash total for the SOURCE PROGRAM deck statements as 
listed in this write-up is 

28381823844211906793 

(2) Hash total for the PROGRAM DECK, consisting of 582 
cards, is 



53313482029556376442 



